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ABSTRACT 

A Study of a Dual Polarization Laser Backscatter System 
For Remote Identification and Measurement 
of Water Pollution. (May 1974) 

Thomas Carlyle Shelves, B.S., Texas A§M University 
and B.A., Baylor University 
Chairman of Advisory Committee; Dr. J. W. Rouse, Jr. 

This study examines the applicability of a dual pola- 
rization laser backscatter system (lidar polarimeter) for 
remote identification and measurement of subsurface water 
turbidity and oil on water. Analytical models for de- 
scribing the backscatter from turbid water and oil on 
turbid water are developed and compared with experimental 
data. Lidar measurements from natural waterways are also 
presented and compared with ground observations of several 
physical water quality parameters. 

The analytical model for describing the backscatter 
from smooth surfaced, turbid water includes subsurface 
single scatter and multiple scatter effects. The like 
polarized backscatter is modeled by an incoherent sum of a 
subsurface single scatter term and a subsurface multiple 
scatter term. The single scatter term is represented by 
a Mie scatter model, incorporating the Mie scattering and 
extinction coefficients. The multiple scatter term employs 



a multiple scatter volume reflection coefficient, similar 
to that used in the analytical treatments of backscatter 
by Rouse [30] and Wilhelmi [18]. Experimental laboratory 
data presented in this report verify the validity of the 
analytical model and demonstrate the characteristic that 
the depolarization ratio (cross return divided by like 
return) increases with an increase in water turbidity. 

Field measurement data from several natural waterways 
show that the like and cross returns, and the depolariza- 
tion ratio, vary directly with the water quality para- 
meters turbidity and suspended solids, and inversely with 
transmittance. These data support the use of a lidar 
polarimeter for remote measurements of these physical 
water quality parameters. 

The analytical model for backscatter from oil on 
water includes the effect on the backscatter of the change 
in refractive index and the attenuation due to oil. The 
oil layer is considered to be a laterally inhomogeneous, 
lossy dielectric. Experimental data performed in the lab- 
oratory verify the validity of the analytical model and 
establish the following characteristics of the data: 1) 

the presence of oil on turbid water does not significantly 
alter the depolarization ratio; and 2) the attenuation due 
to oil is a function of oil type, oil thickness, and exci- 
tation wavelength. Considering these effects on the 
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backscatter, a detection scheme is presented whereby the 
presence of oil on water can be established by using a 
single wavelength lidar polarimeter. A detection, type 
identification, and thickness measurement method is also 
presented using a dual wavelength lidar polarimeter. Using 
these results, the applicability of a lidar polarimeter for 
oil detection, identification, and thickness measurement 
is established. 
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CHAPTER I 
INTRODUCTION 

The water environment can generally be character- 
ized as a dilute, aqueous solution, containing a large 
variety o£ dissolved and suspended, organic and inorganic, 
chemical species, and including an abundance of plant and 
animal life [1]. Because of the invaluable productivity 
and aesthetic value that this environment provides, the 
control of industrial, agricultural, and municipal water 
waste discharges has become increasingly important. Con- 
trol efforts exerted by governmental agencies have been 
segmented into two areas: the control of surface efflu- 

ents, and the control of subsurface effluents. The most 
harmful surface effluent presently receiving waste control 
attention is oil on water. Several ecological investiga- 
tions have shown that petroleum products discharged or 
spilled in the water environment present a serious threat 
to the existence of the aquatic habitat. Suspended sub- 
surface pollutants from precipitation run-off or industri- 
al discharge have also generated significant concern among 



2 


environmentalists because of impairment to the aquatic 
habitat and degradation of the aesthetic value. 

While governmental agencies have attempted to 
control surface and subsurface pollution, their efforts 
have been handicapped by the ineffectiveness of present 
day instrumentation methods for identifying and measuring 
water pollutants. These methods usually require either 
slow laboratory analysis, or expensive field equipment 
analysis. This report discusses some of the present 
methods of water quality analysis, and then presents a 
new method utilizing laser backscatter for performing 
remote, continuous, real-time identification and measure- 
ment of oil on water and selected water quality parameters. 

Water Quality Determination 

The purpose of water quality parameters is to 
provide a means for quantitatively describing water 
impurities. Three basic classifications of water quality 
parameters are biological, chemical, and physical, each of 
which refers to the nature of the water impurity. This 
section describes some of the present methods for deter- 
mining these parameters. 
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Direct Sampling Methods 

The traditional approach to measuring water quality 
parameters involves direct sampling methods. These mea- 
surement methods can generally be divided into two cate- 
gories: 1) laboratory measurements and 2) in situ 

measurements. Laboratory measurements involve collection 
of water samples from the sites, followed by an analysis 
of the samples in the laboratory. There are two major 
problems associated with this approach: 1) the analysis 

of the samples is not real-time and therefore cannot 
adequately monitor the dynamic pollutant variations that 
occur in the water environment; and 2) the physical and 
chemical form of the desired constituent may be altered as 
a result of moving the sample from its natural environment. 

In situ measurements attempt to alleviate some of 
these problems. This measurement method requires portable 
field equipment which can determine water quality para- 
meters at the site. This method is usually preferrable to 
the laboratory method but still has several inherent prob- 
lems : 1) sample collecting and instrumentation operation 

usually require several qualified personnel; 2) large 
areas cannot be measured in short periods of time; 

3) continuous analysis is usually not possible with present 
instrumentation methods; and 4) the measurement procedures 
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are usually very tedious. 

The realization of the need for continuous, real-time 
water quality analysis has motivated several investigators 
to consider remote sensing techniques for water quality 
monitoring. The following section describes some of these 
investigations . 

Remote Sensing Methods 

The basic approach of remote sensing is the detec- 
tion and measurement of electromagnetic radiation for the 
identification and measurement of geophysical phenomena. 

The source from which the electromagnetic radiation origi- 
nates defines two types of remote sensing systems: passive 
systems and active systems. Passive systems identify and 
measure geophysical phenomena through the detection and 
measurement of naturally occurring radiation, such as 
solar or thermal radiation. Active systems, however, 
generate their own source energy and identify or measure 
geophysical phenomena through the measurement of the 
scattered, reflected, or transmitted source electromagnetic 
radiation . 

Passive Systems . The simplest and oldest form of 
passive remote sensing is photography. The primary advan- 
tage of aerial photography is the performance of high 
ground resolution mapping with large areal coverage. For 
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determining water quality parameters, aerial photography 
senses natural radiation scattered or emitted from the 
water subsurface. Several aerial photography techniques 
have been used to detect water pollution, including the 
use of color film, infrared color film, and selective 
spectral filtering [2], 

One problem that has restricted the successful 
application of aerial photography for water quality analy- 
sis has been the inability to separate peripheral atmos- 
pheric effects from the desired water subsurface radiation. 
Ground control reflectances have been used to suppress 
peripheral atmospheric effects but these have provided 
only marginally satisfactory results. Piech, e^ al^. [3,4], 
however, introduced a scene color standard technique which 
provides accurate photographic measurements of volume spec- 
tral reflectances without interference from atmospheric 
effects. Their study showed the feasibility of accurate, 
quantitative analysis of physical water quality parameters 
through aerial photography. 

Other studies in aerial photographic analysis of 
water quality have been performed by James and Burgess [5]. 
Through a computer analysis of data obtained from densi- 
tizing photographic data, an analysis of waste concentra- 
tions of pulp mill outfalls was made with the results 
giving approximately 90% correlations between photographic 
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data and ground observations of several physical water 
quality parameters. 

In 1971, Halajian and Hallock [6] introduced a direct 
digital approach to the imaging and signature analysis of 
water targets and indicated that quantitative analysis of 
water turbidity and water depth was feasible. Later in 
1971, Hallock and Halajian [7] reported a similar investi- 
gation using a passive dual-polarization system for remote 
sensing of water quality. This study indicated that dual 
polarization measurements offer another means of separating 
peripheral atmospheric effects from measurements of water 
subsurface reflectance. 

Schwebel [8] has also shown the feasibility of deter- 
mining turbidity using photographic methods. Color pho- 
tography with selective spectral filtering was taken of 
numerous ponds with simultaneous measurements of ground 
observations of water turbidity. Correlations between 
90% and 95% were obtained from a regression analysis of 
the photographic data with water turbidity. The useful- 
ness of this method, however, was restricted by 1) the 
use of ground control reflectances, and 2) day-time opera- 
tion only. 

Passive infrared imagery has been used by Foster [9] 
to study thermal mixing in effluent waters. The results 
showed good accuracy in quantitative remote measurement 
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of water temperature. Thermal infrared imagery has also 
been applied to oil slick detection [10]. The imagery 
showed that areas with thin layers of oil appeared darker 
[colder) than the surrounding water, while areas with thick 
layers of oil appeared lighter (warmer) . 

Although some passive remote sensing systems have 
performed quantitative analysis of a few water quality 
parameters, these systems have had several disadvantages: 

1) sensing is generally restricted to daytime only, 

2) sensing is generally restricted to clear weather; and 

3) parameter analysis is generally not real-time. 

Active Systems . Active remote sensing is a relatively 
new approach for identifying geophysical characteristics. 

The primary advantage of active sensing over passive sens- 
ing is the capability of: 1) nighttime operation, and 
2) detection insensitivity to ambient radiation. 

Active microwave (radar) systems have been used in 
many earth resource applications. A recent investigation 
in water quality analysis includes the use of a four fre- 
quency microwave radar system for detecting oil slicks [11]. 
The results of this study indicate that detection of oil 
thickness of 1 pm are feasible, although classification 
of oil type could not be established. 

The advent of the laser has provided a powerful new 
tool for application in active remote sensing of water 
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quality. Kim [12] has studied the application of an air- 
borne laser fluorosensor for the detection of algae in 
the sea and has demonstrated the capability for measuring 
chlorophyll-a concentrations in a range from 0,0 to 
30 mgs/m . Laser flourosensors have also been used for the 
detection of oil spills and have shown the feasibility of 
classifying oil types and thicknesses [13]. 

Granastein [14] reported in 1972 a laboratory study 
which used a linearly polarized laser and a beaker of water 
with known suspended scatterers and absorbers to determine 
the relation between the volume scatter and the concentra- 
tions scatterers and absorbers. Polarization effects were 
also considered in a later study [15], The results indi- 
cated that the depolarization ratio (orthogonally depol- 
arized backscatter divided by like polarized backscatter) 
depended on the ratio of the scattering to absorbing con- 
centrations, and not the absolute mass concentrations. 

This work is significant in that it shows the dependence 
of the depolarization on subsurface effects. 

Another study investigating laser backscatter from 
turbid water was conducted by Wilhelmi, Mayo, and Rouse 
[16]. A dual polarization laser backscatter system (lidar 
polarimeter) was constructed and used to measure backscat- 
ter from laboratory simulations of natural water. The 
hypothesis which motivated the design and construction 
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of the lidar polarimeter was that suspended contaminants 
in real water scatter and depolarize the incident radia- 
tion, and that the degree of depolarization should indicate 
concentrations of contaminants. Laboratory measurements 
of backscatter supported this hypothesis, and therefore 
the investigation established the feasibility of deter- 
mining certain water quality parameters with a lidar polar- 
imeter . 

A further investigation with the lidar polarimeter 
was performed by Shelves [17] . In this study, an analyti- 
cal model was presented which described the backscatter 
from smooth surface turbid water illuminated by vertically 
polarized laser light. Experimental measurements were 
performed and compared with calculated values. The results 
of this study were: 1) a further indication of the feasi- 

bility of using a lidar polarimeter for determining water 
quality parameters, and 2) the postulation that a volume 
multiple scatter mechanism should explain the depolariza- 
tion process. 

The most recent investigation with the lidar polari- 
meter was performed by Wilhelmi [18]. This study included 
a theoretical development of the polarized and depolarized 
backscatter component from a rough surface and a smooth 
surface, incorporating single scatter surface and multiple 
scatter subsurface effects. Experimental laboratory 
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measurements were performed to investigate the validity of 
the theoretical model. The results of this study indicated 
that since the theoretical values and the experimental 
values had correlation coefficients greater than 90%, a 
volume multiple scatter mechanism appears to explain the 
depolarization process. 

Report Objectives 

The foregoing discussion shows two important results 
which provide the stimulus for this investigation. First, 
it establishes the need for real-time identification and 
measurement of surface and subsurface water quality para- 
meters; and second, it demonstrates the feasibility of 
using laser backscatter for this application. This inves- 
tigation extends the previous laser backscatter studies 
in the following manner. First, a new analytical approach 
is taken to describe the backscatter from turbid water. 
Previous laser backscatter studies, such as that performed 
by Wilhelmi [18], assumed that the backscattered energy 
was entirely multiple scatter from the subsurface. 

However, backscatter measurements from turbid water show 
that the like-polarized component is significantly larger 
than the cross-polarized component which strongly suggests 
the presence of single scatter within the volume. 

Therefore, the model presented in this report incorporates 
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a single scatter component into the backscattered return. 
The results of laboratory measurements are also given to 
show the validity of the model. 

A second extension of previous studies presented in 
this report is the development of an analytical model for 
describing the backscatter from oil on turbid water. An 
identification and measurement method using a lidar polari- 
meter for detecting oil on water is discussed along with 
preliminary laboratory measurements to examine the concept. 

Furthermore, this report describes the results of 
numerous lidar polarimeter field measurements conducted 
from several natural waterways. A discussion of the field 
portable lidar polarimeter used in these field measurements 
along with an analysis of the data is included. 

Scope of Report 

Chapter II provides the necessary background for 
understanding the nature of water pollution and for selec- 
ting water quality parameters which appear suitable for 
detection by the lidar system. 

Chapter III describes the analytical modeling tech- 
niques used to characterize the electromagnetic interaction 
between the water subsurface and the lidar system. Also 
included is an electromagnetic scattering model for oil 


on water. 
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Chapter IV describes the lidar polarimeter used in 
this investigation including a description of the system 
components and transceiver geometrical configuration. 

Chapter V describes the experiments of this study. 
Lidar laboratory measurements with simulated turbid water 
and oil on water are discussed in conjunction with lidar 
field measurements performed from several natural waterways. 

Chapter VI compares the experimental results with 
the analytical models for backscatter from turbid water 
and oil on turbid water. Lidar field measurement data 
are also presented graphically to show the relationships 
between the lidar returns and several water quality para- 
meters . 

Chapter VII gives the conclusions drawn from this 
study concerning the validity of the analytical models and 
the applicability of using laser backscatter for deter- 
mining water quality parameters. In addition, this 
chapter provides recommendations for future investigations. 
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CHAPTER II 
WATER QUALITY 

This chapter discusses some of the characteristics 
of the water environment and concludes with a selection 
of surface and subsurface water quality parameters which 
appear applicable to detection and measurement by the 
laser sensor. 


Water Environment 

The water environment is a very complex, dynamic, 
heterogeneous system. The structure of the water environ- 
ment generally consists of four constantly exchanging 
aquatic systems: estuaries, rivers, lakes, and oceans. 

The most dynamic and complex aquatic system is the estuary, 
defined as the mixing basin where the ocean water and 
fresh water rivers meet. The estuary is characterized by 
extremes in water quality, currents, and bottom deposits 
[19]. In addition, the intertidal zone is subjected to 
alternate exposure of air and water. 

The fresh water rivers and lakes are burdened with 
the transport of waste discharges and sediment deposits. 

The ability of the rivers to assimilate these wastes 
depends largely on the type and volume of the wastes. 

When the levels of waste discharges are low, the water 



removes the wastes and cleans itself. At higher levels of 
waste discharges, the ability of the rivers to assimilate 
wastes becomes restricted, resulting in environmental 
degradation . 

The oceans, however, have almost an unlimited ability 
to assimilate wastes, and because of their relatively con- 
stant characteristics, absorb most of the effluents 
discharged into the estuary by the fresh water rivers. 

The dynamic nature of the aquatic environment is a 
result of several factors, besides the obvious factor of 
varying waste discharges. Climatic conditions such as 
runoff from snowmelt or floods can produce muddy, soft 
water with a high bacterial count. Runoff during a 
drought, however, can produce high mineral content 
groundwater . 

Geographic conditions also are a primary factor in 
water quality variations. For example, in the great belt 
of heavily populated areas along the east side of the 
Mississippi, sediment concentrations in most surface 
waters average less than 270 mg/liter. For the eastern 
United States though, sediment concentrations average 
approximately 1900 mg/liter [20]. 

Seasonal conditions such as aquatic organism growth, 
floods, droughts, waste discharges, and the overturn lakes 
and reservoirs, also contribute to the dynamic character 
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o£ the water environment. 

Physical Subsurface Impurities 
General Characteristics 

Subsurface impurities in the aquatic environment, 
whether solid, liquid, or gas, are dispersed into three 
progressively finer states: suspended, colloidal, and 

dissolved. Dissolved impurities are those which are dis- 
persed in water as single molecules or ions [20]. Gener- 
ally, particles less than .01 ym are classified as 
dissolved impurities. 

Insoluble material is usually defined as any particle 
from .01 to 200 ym in diameter, and includes both suspended 
and colloidal materials. Particles greater than 200 ym 
will not generally be found suspended in water unless the 
water has high velocity currents, or unless the matter is 
bouyant bio-organic particles [20] . 

Many investigations have been conducted to study the 
size distribution of particles suspended in the ocean. 
Sheldon, et al. [21] used the Coulter counter technique to 
study suspended particles in the North Atlantic Ocean sur- 
face waters and found particle sizes with geometric mean 
diameters of approximately 20 ym. Gordon [22] also con- 
ducted a similar type study in the North Atlantic Ocean 
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and found that particles with diameters less than 7 pm 
were the most numerically abundant and that there was no 
significant change in size distribution with ocean depth. 

Krey [23] conducted a study in estuaries and found 
the mean spherical diameter of suspended organic particles 
to be between 5-6 pm. Biggs conducted a sediment size 
analyzer of the Chesapeake Bay and reported sediment size 
distributions with geometric means of approximately 25- 
30 pm. The significance of these example particle size 
distributions will become apparent when analytical modeling 
^ is considered in Chapter III. 

Measurement Parameters 

Physical water quality parameters refer to the 
appearance, taste, or odor of the water and thus are a good 
measure of the aesthetic quality of the water. The appear- 
ance of the water is usually measured by light scattering 
techniques. Most of these techniques consist of a colli- 
mated beam of light illuminating a water sample with the 
scattered, or transmitted light collected by a light 
detector which may be spectrally filtered. The intensity 
of scattered or transmitted light is compared with a known 
standard, with the results giving a quantitative measure 
of the appearance of the water. 
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Several physical subsurface water quality parameters 
appear suitable for measurement by the laser sensor 
examined in this report. These parameters include turbid- 
ity i suspended solids, and transmittance, and are discussed 
in the following paragraphs. 

Turbidity . Turbidity is a measure of the light 
scattering characteristics of water caused by the presence 
of colloidal and suspended matter. The presence of sus- 
pended matter can indicate a change in the water quality, 
for example, contamination of the water sample by micro- 
organisms or the presence of finely divided inorganic 
substances such as clay or silt. The determination of 
turbidity affords a very sensitive means of detecting the 
presence of suspended matter. Of some importance is the 
fact that the presence of suspended matter will degrade 
the appearance of the water. Turbidity is also undesira- 
ble for industry applications where the product is destined 
for human consumption, for domestic water supplies, and 
for other industrial applications such as pulp and paper 
manufacture. The turbidity of natural waters is also an 
important factor in the control of productivity. Turbidity 
interferes with the penetration of light and will affect 
the ecosystem since photosynthetic activity within the 
ecosystem is highly dependent upon light penetration. 
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The traditional procedure for the measurement of 
turbidity is the visual candle turbidimeter, usually 
termed the Jackson candle turbidimeter. In this method 
the extinction of the image of a candle flame, situated 
below a vertically suspended graduated tube, is observed 
by viewing the sample through various depths of the water. 
The sample tube is graduated in Jackson Turbidity Units 
(JTU) . The lower limit of turbidity that may be measured 
using the Jackson candle turbidimeter is 25 JTU. 

A new measurement technique for measuring turbidity 
termed nephelometry has been recommended by Standard 
Methods [25]. A nephelometer measures the intensity of 
scattered light at right angles to the incident collimated 
light beam. The measurement units of a nephelometer are 
Formazin Turbidity Units (FTU, where one FTU is approxi- 
mately equal to one JTU) named after the standard solution 
which is used to calibrate the meter. The primary advan- 
tages of nephelometers are: 1) the ability to measure 

turbidities much less than 25 FTU; and 2) the elimination 
of visual judgments as required by the Jackson candle 
turbidimeter. 

The sampling frequency for measuring the turbidity 
of an estuary or stream depends on how critical the 
quality of the water is to the particular location. How- 
ever, as a general rule, if turbidity is an important 
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parameter, sampling should be continuous [1] . 

Suspended Solids . Suspended solids is a measure of 
the weight concentration of particulate matter in water. 
Although the presence of suspended solids does cause tur- 
bidity, the amount of turbidity in the water does not 
necessarily relate to the amount of suspended solids. The 
measurement of suspended solids is made by filtering a 
known volume of water and determining the weight difference 
between the clean filter and the contaminated filter. The 
units of suspended solids are usually mg/liter. 

Transmittance . Transmittance is a measure of the 
light attenuation of water, and is related to the scatter- 
ing and absorbing constituents in the water. Transmittance 
is usually measured with a spectrally filtered light source 
and detector, with units in percent, relative to distilled 
water. The importance of transmittance in ecology is 
similar to that of turbidity. 

Oil on Water 

Oil pollution has become a problem of major impact 
on the environment - a problem of a more pervasive and 
disturbing nature than the obvious detrimental effects 
observed on the habitat and shores. Before examining these 
detrimental effects, some of the general characteristics 
of oil pollution are discussed. 
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General Characteristics 

Oil is not a single substance but a complicated and 
variable mixture of hundreds of chemical compounds. The 
constituents of oil share many common properties but also 
differ considerably in their influence on the environment 
[26] . Among these properties are: 

1. Toxicity - Many low boiling aromatic hydrocarbons 
are lethal poisons to almost any organism, while 
some higher boiling paraffin hydrocarbons are 
essentially nontoxic to most forms of life. No. 2 
fuel oil, for example, is widely known to be 
extremely lethal to marine life while in general 
crude oils are not highly toxic [27]. 

2. Solubility - Petroleum derivatives in certain 
concentrations are soluble, while other hydro- 
carbons are essentially insoluble. Solubility 
significantly influences the toxicity of a 
component of oil. 

3. Biodegradability - This varies widely according 
to the nature of the hydrocarbon with the rate 
of biodegradability significantly effecting the 
persistance of environmental effects. 

4. Carcinogenity - Some components of oil are known 
to have cancer- inducing properties [26]. 
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The effect of oil on the environment' is not only 
dependent on the type of oil spilled as discussed above, 
but also on the amount spilled. Typical amounts of oil 
spills range from a few hundred barrels to hundreds of 
thousands of barrels, as occurred in the famous Torrey 
Canyon spill where 400,000 barrels of oil were spilled. 

The spreading of oil from such spills largely depends 
on the type of oil. Pollution by light oils in the 
estuaries (where unfortunately most spills occur) seldom 
reach the shore. Rather, they dissipate quickly, parti- 
cularly if disturbed by wind or waves. Small scale spills 
of crude oils, however, tend to move inland. In large 
scale crude oil spills, the oil remains in a cohesive 
mass and comes ashore in the form of a slick. Pollution 
by heavy oils tends to breakup into large lumps and, parti- 
cularly in cold weather, solidifies very quickly [28]. 

Thicknesses of oil spills after spreading vary 
largely with the condition of the wind and waves. Gener- 
ally, however, thicknesses for most spills range from 
tenths of millimeters to several millimeters. 

Effect on Aquatic System 

The overall effect of oil pollution on the aquatic 
system is best illustrated by the following oil spill 
example discussed by Murphy [26]. The example discussed 
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is the West Falmouth spill which occurred near the shores 
of Cape Cod in 1969 with 4,000 barrels of No. 2 fuel oil 
spilled. A feature of this medium size spill was a massive 
initial kill of area marine life, in accordance with the 
nature of most spills of this type. Most of the spill 
dispersed into the sea after two days, resulting in a 
persistence of oil on the ocean floor which expanded to 
cover 50 hectacres (more than 1 acre per original barrel 
of oil spilled). As the spill expanded on the ocean floor, 
the mortality to bottom life increased significantly. The 
kill of marsh plants and bottom organisms apparently 
changed the physical properties of the sediments and may 
have caused erosion and spreading of the trapped oil. 
Another feature of this oil spill was the tainting of the 
flesh of commercial fish and shellfish, which has been one 
of the oldest and most frequent complaints of oil pollu- 
tion. 

The effects of oil spills on the aquatic system vary 
for different oil spills. However, the West Falmouth 
spill is described as a typical oil spill and is considered 
to be representative of the impact of oil pollution on the 
aquatic system. 
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CHAPTER III 

ANALYTICAL MODEL FOR LASER BACKSCATTER 

Many investigations have been conducted to determine 
an analytical model for describing backscatter from a 
volume of scatterers. Some investigations have presented 
single scatter Mie models to describe the backscatter, 
however, numerous laboratory measurements with spherical 
scatterers illuminated by a polarized light source have 
shown a significant amount of depolarization of the back- 
scattered energy, which suggests the presence of multiple 
scatter within the volume. While several subsurface mul- 
tiple scatter models have been investigated for describing 
backscatter from a volume of scatterers [29] , [30] , these 
models assume that the backscattered energy is randomly 
polarized, that is, the backscattered electric field has 
no preferred polarization (For a comparison of these 
scattering models, see Rouse [31]). Measurements from 
cast dielectric samples and turbid water reported by 
Wilhelmi [18], however, show large like polarized (rela- 
tive to cross polarized) returns, which suggests the 
presence of single scatter within the volume. This chapter 
presents the development of a volumetric backscatter model 
which incorporates both a single scatter return and a mul- 
tiple scatter return from an aqueous suspension of 
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scatterers illuminated by a polarized electromagnetic 
plane wave. These results are then extended to include 
the effect on the backscatter of a layer of oil on the 
turbid water. 


Turbid Water 


Approach 

The following approach will be taken in developing 
the volumetric backscatter model. A linearly polarized 
electromagnetic plane wave is incident upon a smooth sur- 
face volume containing spherical scatterers (Figure III-l). 
Since the backscatter from a smooth surface can be assumed 
negligible for optical frequencies [18] » these effects 
will not be considered. The backscatter from the subsur- 
face will consist of two components: a like polarized 

return and a cross depolarized return. The like polarized 
return (P^y oi" Pj^j^* where the first subscript indicates 
transmitter polarization, vertical or horizontal, and the 
second subscript indicates receiver polarization) consists 
of a single scatter return (denoted by superscript s) plus 
a multiple scatter return (denoted by superscript m) . 

Thus , 
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Figure III-2. Simplified Scattering Geometry 
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The cross depolarized return or from spherical 

scatterers consists only of a multiple scatter return or 




(III-3) 




(III-4) 


The following sections develop this model for single 
scatter return and the multiple scatter return. 

Single Scatter Model 

Assuming a vertically polarized, monochromatic plane 
wave incident upon a smooth surface volume of water and 
considering the simplified geometry shown in Figure III- 2 
(since the incidence angle 6 shown in Figure III-l effects 
only the transmission coefficients) , the irradiance at 

the top surface of the water is 


Hys = /'^s 



27 


where is the power transmitted (assuming lossless trans- 
mission) and is the area of illumination at the surface 
s 

of the water. The portion of the irradiance trans- 

^ vw 

mitted into the water is 

Hyw 

where is the Fresnel power transmission coefficient 
for vertical polarization for propagation from medium 1 
(air) to medium 3 (water). 

Due to the attenuation of the beam through the volume 
(due to scattering and absorption by the scatterers) , the 
irradiance H ^ on a differential element within the volume 

VO 

is 

Hvo= Hyyy, exp(-^x) 



where g is the volume extinction coefficient and x is the 
distance from the surface to the differential volume. The 
volume extinction coefficient is defined by van de Hulst 
[32] as 


^ ■ / N(a)-7ra^da 


(in-5) 


where is the efficiency factor for scattering and 

absorption (polarization independent for spherical parti- 
cles), and N(a)da is the number of particles per unit 
volume with radii a in the interval da. For a suspension 
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of particles with a broad size distribution (polydisperse) , 
it is advantageous to express the volume extinction coeffi- 
cient in terms of the total number of particles per unit 
volume by defining 

N(a)da = Nr(a)da 

where N is the total number of particles per unit volume 
and r(a)da is the relative number of particles per unit 
volume with radii a in the interval da. Equation (III- 5) 
can then be written as 

^ = nJ Qgxt(a)'^a^i'(a)da 

Using this equation to calculate 6 requires only one 
numerical evaluation of the integral expression for differ- 
ent particle concentrations. Defining this integral to be 
3^ gives 


^ (III-6) 

For a suspension of particles with a very narrow size 
distribution (monodisperse) , the volume extinction coeffi- 
cient is simply 


^ ■ ^ext N TT 
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The general relationship for the irradiance scattered 
from a volume dV in the direction is 

dH = H,F(*)dV/ 

where is the scattering function, which in general is 

related to the polarization state, relative particle size 
distribution, particle orientation, and particle concen- 
tration [32]; ip is the scattering angle measured from the 
direction of the incident wave; r is the distance from the 
differential volume to the point of observation; and is 
the incident irradiance. The scattering function for back- 
scatter from spheres is the same for each of the 

plane polarized components of incident radiations, so that 
it is not necessary to designate the polarization of the 
scattering function [33]. If the relative size distribu- 
tion is used to calculate the scattering function F(7r), 
then F(7 t] can be given by 

FCtt) = N (tt) (III- 7) 

where F^Ctt) is the scattering function calculated from the 
relative particle size distribution. 

Assuming that the distance R >> x (Figure III-2, p.25) 
as is the case for most applications, the distance r can be 
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approximated by R. The irradiance backscattered from a 
volume dV is then 

dHyy=l^(7]g Tg'^) NF^(ir)exp(-2N^fX)dV cm-s) 
AsR2 

where is the Fresnel transmission coefficient for 
vertical polarization for propagation from medium 3 to 
medium 1. Since single scattering is assumed, the volume 
element dV must be within the boundaries of the incident 
beam and can thus be represented as 

dV=A^dx 

o 

Substituting this equation into (III- 8) and confining 
the limit of integration on x to be the depth L along the 
beam gives 

NF^W / exp(-2No,x) dx 

0 . r2 0 

which after integrating gives 


^vv" Tsi) W ( I - exp{- 2N L )) 

r 2 Zlv 
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The single scattered power collected by a 

receiver o£ area is then 

r 

= Ta^) l?(7r)(l-exp(-2N^rL)) (m.9, 

r 2 2^r 

where the factor n is a system constant to correct for the 
efficiency of the optical components and the calibration 
of the system. 

It is appropriate at this point, to define for the 
single scatter model an effective volume reflection coeffi- 
cient Observing (III-9), the terms which describe the 

scattering media in terms of its scattering properties 
are the scattering function F^(tt) and the volume extinction 
coefficient 3^. The combination of the scattering func- 
tion and the volume extinction coefficient derived through 
the integration process represents the fraction of the 
incident energy which is single scattered to the receiver. 
This fraction of energy will be defined to be the single 
scatter volume reflection coefficient which for verti- 
cal transmit polarization is given by 

\> l?(-n-)(l-exp(-2N/3,L)) 


Several observations can be made about (III- 10) con- 
cerning the magnitude of the volume reflection coefficient 



relative to the mean depth L and the volume extinction 
coefficient. Representing the exponential in CIH"10) 
in a series expansion gives 
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exp(-2N = I - 2N + (2N^Lf ♦ 

2! 

Assuming that N3^L << 1 (or that the mean depth is much 
smaller than the extinction length, 1/N3^) gives 

exp(-2N^f L) I -2N/8^L 

Substituting this expression into (III- 10) gives 

(iii-ii) 

which shows that at very low particle concentrations, the 
single scatter return is proportional to the particle con- 
centration (consistent with single scatter theory) , 

Considering the case where Ng^L >> 3 (or that the 
extinction length is 3 times smaller than the mean depth 
L) , the single scatter volume reflection coefficient is 
approximated by 


(III-12) 
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which shows that the single scatter return is independent 
of particle concentration. 

For the case of Np^L==l, (III- 10) must be used. 

Multiple Scatter Model 

Several assumptions concerning the nature of the 
multiple scatter mechanism will be made in considering an 
analytical model for multiple scatter within the volume. 
First, no attempt will be made to provide a physical model 
of the multiple scatter mechanism. Instead, the more 
general approach is used by assuming that the backscatter 
from a multiple scattered wave exists which is some frac- 
tion of the incident field. Furthermore, the scattering 
media is assumed to be bounded by a smooth surface. 

Following the same procedure for the multiple scatter 
case as for single scatter and noting that the multiple 
scatter term should be of the same form as (III-9) , a 
volume reflection coefficient is used to represent the 

fraction of the backscattered energy which is like- 
polarized multiple scatter. The backscattered power for 
this case can then be given as 



r2 


(III-13) 
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Realizing that due to multiple scatter there is a 
significant amount of energy contained within the cross 
polarized component, the same procedure is followed by 
defining a volume reflection coefficient The volume 

reflection coefficient is treated as a tensor to account 
for any polarization dependency of the particle scattering 
within the media. The cross polarized power can then be 
given as 

Pyh =1? Ar(T,^T3*?) c:ii-i4) 

Note that the transmission coefficient for the cross term 
contains the term T^j^, which accounts for the cross polar- 
ized transmission from media 3 to 1. 

An important observation concerning the previous 
developments is that the effective volume reflection coef- 
ficient for the single scatter term is given in terras of 
the physical scattering properties of the medium, that is, 
in terms of the scattering function and extinction coeffi- 
cient. The volume reflection coefficient for the multiple 
scatter process, however, is not given in terms of the 
media properties and must then be determined empirically 
or analytically by another approach. 
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Results 

The model describing the backscatter from a suspension 
of spherical particles bounded by a smooth surface and 
illuminated by a vertically polarized plane wave, can be 
obtained in terms of the like polarized components and 
cross polarized components as 


R2 



CIII-15) 



(III-16) 


Similarly, for horizontal transmit polarization the like 
and cross components are 


hT-h. 


R* 


CIII-18) 
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Oil on Water 


Approach 

Oil on turbid water can be described as a laterally 
inhomogeneous, isotropic, lossy dielectric. (It is assumed 
that the oil is homogeneous with depth, that is, the index 
of refraction of the oil is not a function of oil depth.) 
The effect of oil on water is twofold: 1) the difference 

in the index of refraction of oil relative to water alters 
the transmission coefficients at the air/oil boundary; and 
2) the lossy properties of the oil results in attenuation 
of the energy transmitted through the oil. The following 
sections extend the turbid water model developed in the 
previous section to include these two effects of oil on 
water. 

Transmission Coefficients 

The effect of a layer of oil on water is diagrammed 
in Figure III-3. Medium 2 can be described as a dielec- 
tric mixture of oil and water since the oil is not later- 
ally homogeneous. To consider this inhomogeneity, an 
effective index of refraction of the oil and water mixture 
is used by weighting each index of refraction by the frac- 
tional areal coverage of oil. Defining this fraction of 
areal coverage to be p, the effective index of refraction 
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n of the oil and water mixture can be given in terms of 
ow 

the relative indices of refraction of oil n^ and water n 

0 w 

as 

® ^ % * (l-p)n^ ciii- 19 ) 

To account for the transmission of energy across boundaries 
1-2 and 2-3 requires the use of two transmission coeffi- 
cients. Two more transmission coefficients are then 
required to account for the transmission of energy out of 
the water volume and across boundary 2-1. 

Attenuation 

The second and more profound effect on the backscatter 
from oil on water is the two way attenuation characteris- 
tic which reduces the backscattered power exponentially 
with respect to the oil thickness, and the oil extinction 
coefficient. Expressing this effect analytically and 
including the boundary transmission effect, the results of 
the turbid water model are modified by 

~ (TjgTggTggTgi ) exp(- 2 °t)F^y (111-20} 

and 
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Pvh ■ IT,$T,»4T2';)exp(-2.t)P,^ 

(III-21) 

where and are the subsurface backscattered powers 
in the absence of oil, given in (III-15) and (III-16); the 
T.^'s are the Fresnel transmission coefficients where the 

ij 

superscript x represents the polarization and the sub- 
scripts represent from medium i to medium j; a is the 
effective extinction coefficient (function of oil type 
and wavelength); and t is the oil thickness. The trans- 
mission coefficients in the denominator normalize P^^ and 
to cancel the air to water transmission coefficients 
observed in (III-15) and (III-16). The effective extinc- 
tion coefficient can be expressed in terms of the extinc- 
tion coefficient for a homogeneous oil layer by the 
fractional cover p as 

a = a.p (III-22) 

n 

Since the oil on water model includes the use of an 
effective refractive index and extinction coefficient, 
one primary restriction of the applicability of model is 
that a long term time average of the return signal and/or 
a large spatial average of the surface must be used so 


' 'i 3 '31' 
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that the effective coefficients mentioned above are appli- 
cable to the measurement situation. 
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CHAPTER IV 
LIDAR POLARIMETER 

The sensor used in this investigation was a dual 
polarization laser backscatter system termed a lidar 
polarimeter. This chapter describes the system components 
[16] and the transceiver configuration of the lidar 
polarimeter. 

System Description 

The lidar polarimeter transmits a linearly polarized 
laser beam and receives both the vertical and horizontal 
polarized backscattered energy. This section describes 
the system components which perform this task. 

Transceiver 

The transmitter of the lidar polarimeter (See 
Figures IV- 1 and IV- 2) is a 5 mw, helium- neon, cw laser 
(Spectra Physics Model 120) provided with a Spectra 
Physics broadband polarization rotator attached to the 
laser output. The polarization rotator allows variation 
of the transmitter polarization. The transmitted laser 
beam is modulated by a Princeton Applied Research (PAR) 
rotary, light chopper to allow detection insensitivity 
to ambient radiation (See Detection Electronics) . 
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Figure IV-1. Lidar Polarimeter Transceiver 
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The receiver of the lidar polariraeter includes 
several separate standard photographic extension tubes 
mated with a Vivitar 135 mm f/3.5 telephoto camera lens. 

The extension tubes house several optical components. A 
Jodon pinhole assembly located in the focal plane of the 
collecting lens provides a narrow receiver field of view, 
with fine adjustment of direction. A five centimeter 
focal length lens recollimates the spatially filtered 
return before passing through a narrow band spectral 
filter. A Spectra Physics 2.5 centimeter beam-splitting 
polarizer cube is located inside an aluminum box assembly 
with ports threaded to fit the extension tubes. The hori- 
zontal and vertical polarization components are detected 
by separate United Detector Model 500 photodiode/op-amp 
combination detectors mounted in extension tubes. The 
characteristics of the transceiver section are summarized 
in Table IV-1, 

The transceiver was mounted on a 19*' x 34” x 1" 
aluminum plate. To allow field portability, the plate 
was then mounted to a heavy duty Majestic instrument 
tripod (Model T-114) with a rotatable head (Model H120) . 

The rotatable head allows continuous tilt from straight 
down to 180“ up. A Davis Sanford tripod dolly Model DS- 
69 was used to aid in sensor mobility. A metal hood 
enclosing the transceiver section assures protection during 



Table IV-1. Lidar Polarimeter Specifications 


Power Transmitted 

Minimum Detectable Polarization 
Ratio 

Receiver Aperture 
Spectral Filter Bandwidth 
Beam Divergence 
Voltage Responsivity 
RMS Noise Deviation 

(1 sec integration time) 
(30 sec integration time) 


6 mw @ 633 nm 
1:200 

38 mm 

3 nm @ 633 nm 
1.7 milliradians 
5 yv/pw 

10 yv 
2 yv 
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transit and inclement weather. Figure IV-3 shows a pho- 
tograph of the transceiver package. 

Detection Electronics 

The modulated signals from the photodiodes are 
synchronously detected with a PAR Model 124 lockin ampli- 
fier, to avoid DC drift, low frequency noise, and ambient 
radiation interference. Measurements with the lockin 
amplifier showed an rms noise deviation of approximately 
10 yv with a 1 second integration time, and 2 yv with a 
30 second integration time. 

The lockin amplifier was housed in a central control 
unit along with the laser power supply and photodiode 
power supply. A description and diagram of the lidar 
polarimeter electronics are given in Appendix A. Figure 
IV- 4 shows the complete lidar polarimeter system. 

Transceiver Configuration 

The bistatic nature of the lidar polarimeter requires 
special considerations. Since the receiver is not on axis 
with the transmitter, an intersection region is created 
at the point where the field of view and transmitted beam 
intersect. Figure IV-5 diagrams this effect. 

The depth of the intersection region can be esti- 
mated by the following procedure. The angle <t> between 




Figure IV-3. Lidar Polarimeter Transceiver 

Package 


48 



Figure IV-4. Complete Lidar Polarimeter 

System 
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Figure IV- 5. Diagram o£ Bistatic Lidar Measurement 
Volume 
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the receiver and transmitter is approximated by 


</> 


W 


R 


assuming that R >> w, where w is the separation distance 
between the transmitter and the receiver, and R is the 
slant range to the surface of the target. Observing the 
geometry of the intersection region and assuming that 
R >> d£ shows that the mean depth L is given by 


or 



(IV-1) 


where d£ represents the diameter of the field of view at 
the surface of the target. Using the ray optics approach, 
the diameter of the field of view at the surface can be 
given in terms of the diameter of the receiver spatial 
filter (pinhole) d^j^ by 



^ph R 
F 


where F is the focal length of the collecting lens. 
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Substituting this expression into (IV-l) gives 


L = 


2 w F 


(IV- 2) 


Another system parameter which will be needed in 
later developments is the diameter of the laser beam dj^ 
at the surface of the target. For the laser used in the 
experiments, dj^ is approximated in terms of the range in 
meters and the beam divergence (1.7 milliradians) as 


cIl - .001 + .00I7R 


The area of illumination can then be given by 


Al= -f-CoOl + ,0017R)^ 


(IV-3) 
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CHAPTER V 
EXPERIMENTS 

Several laboratory turbid water experiments were 
conducted with the lidar polarimeter to determine the 
validity of the scattering model developed in Chapter III 
and to better understand the nature of the scattering 
mechanism. Several field experiments along natural water- 
ways were conducted with the lidar polarimeter to deter- 
mine the applicability of the laser sensor for determining 
subsurface water quality parameters. Laboratory oil on 
water experiments were also conducted to determine the 
validity of the oil on water scattering model and to 
demonstrate the scattering phenomenas involved so that a 
procedure could be presented for characterizing oil pollu- 
tion with respect to oil type and oil thickness. This 
chapter describes these experiments. 

Turbid Water 

Achieving the objectives of the laboratory turbid 
water experiments required varying several parameters. 

These include: 1) varying the suspended particle concentra- 
tion from very low concentrations where single scatter 
should be dominant to very high concentrations where mul- 
tiple scatter should be dominant; 2) varying the receiver 



field of view to determine the dependency of the depolari- 
zation on the scattering volume viewed by the receiver; 

3) varying the transmitter polarization; and 4) varying 
the absorption losses of the water medium to determine the 
effect on the depolarization. 

The first experiments included measuring and 

and varying 1) the suspended particle concentration 
(number of particles, N',/mL) from 1.75 X 10^ to 2.03 X 
10 , and 2) the diameter of the receiver field of view 

at the surface from .32 cm to 1.07 cm. The suspended 
particles used for these measurements were relatively mono^ 
dispersive polystyrene latex spheres (Dow Plastic Pigment 
722) with a known refractive index (1.59) and particle 
size distribution (Figures V-1 and V-2) [34]. The field 
of view was varied by changing the spatial filter diameter 
from .3 mm to 1 mm. A second set of backscatter measure- 
ments with a constant field of view (.64 cm) were made 
by varying the suspended particle concentration using a 
relatively polydispersive polytetraf luoroethylene resin 
dispersion ("Teflon” TFE 30 dispersion manufactured by 
E, I. Dupont Nemours and Co.). This dispersion consisted 
of spherical particles with a refractive index of 1.37 and 
diameters in the range .04 pm to .4 pm (Figures V-3 and 
V-4) [35]. This second measurement also included measuring 
not only P^^ and P^j^ but also Pj^j^ and 
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Figure V-1. Approximate Relative Size Distribu- 
tion for Dow Plastic Pigment 722 



Figure V-2. Photomicrograph of Dow Plastic 
Pigment 722 
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Figure V-3. Approximate Relative Size Distri- 
bution for Dupont TFE 30 



Figure V-4 


Photomicrograph of Dupont TFE 30 


S6 


The last turbid water experiment included measuring 

P and P u from a constant scattering concentration 
vv vh 

(using polystyrene latex particles with a concentration 
of 4.37 X 10^ N’/inL) and a constant field of view (.64 cm) 
with a varying absorbent concentration of the water 
medium. The absorbent used for this experiment was black 
nigrosine dye. Dye concentrations ranged from 0.0 to 
214.0 mg/L. All data collected in these experiments along 
with values for the system parameters are listed in 
Appendix B. 


Natural Waterways 

Several field experiments from natural waterways 
were conducted with the lidar polarimeter from atop 
bridges and boats. The first significant field experiment 
was conducted from a bridge on the Brazos River at Waco, 
Texas (Figure V-5). Lidar measurements of P^^ and P^j^ 
(also a few measurements of P^j^ and Pj^^) were recorded 
over a two day period at an incidence angle of 30'’, a 
slant range of 9.5 m and a reciever field of view at the 
surface of 4.8 cm. Water samples were collected while 
the lidar measurements were being performed so that an 
analysis of turbidity, suspended solids and transmittance 
could be performed for each sample. Special care was 
taken to inhibit physical or chemical alteration of the 
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Figure V-5. Experimental Set Up On 
Brazos River Bridge 


water impurities, A twenty minute time history of the 
like polarized return was recorded with a Hewlett-Packard 
7004-B X-Y Recorder, 

The second field experiment was conducted from the 
boat Excellence on an eight mile stretch along the Houston 
Ship Channel (Figure V-6) . Lidar measurements of and 
made at an incidence angle of 35°, a slant range 
of 2.84 m and a receiver field of view at the surface of 
approximately 1.25 cm. Ground observations of turbidity 
and transmittance were also recorded - 

The final field measurement was conducted again on 
the boat Excellence on a twenty mile stretch along three 
interconnecting waterways: the Intercoastal Waterway, 

the Freeport Ship Channel, and the Brazos River. Lidar 
measurements of P^^ and were made at an incidence 
angle of 35°, a slant range of 2.84 m, and a receiver 
field of view at the surface of approximately 1.25 cm. 
Ground observations of turbidity, suspended solids, and 
transmittance were performed. Data for these field mea- 
surements are given in Appendix B. 

Oil on Water 

The objectives of the oil on water experiments 
included demonstrating four concepts: Ij oil on water 

attenuates the energy backscattered to the receiver; 2) the 
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Figure V-6. 


Experimental Set Up Aboard 
The Boat Excellence 
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presence of oil on water does not significantly alter the 
subsurface depolarization ratio; 3) the extinction coef- 
ficient of a specific oil type is a function of the wave- 
length; and 4) the extinction coefficient varies for 
different oil types. 

The first oil on water experiment conducted with 
the lidar polarimeter included measuring and and 
varying the oil thickness for different test oils. The 
test oils used for these experiments were . gasoline , kero- 
sene, SAE 30 refined motor oil, no. 2 fuel oil, and crude 
oil. Thicknesses ranged from .05 mm to 5.5 mm. The 
subsurface scattering medium for all test oils except 
crude oil was a suspension of Dow Plastic Pigment 722 
particles with a concentration of 1.75 X 10^ N'/mL. A 
much higher concentration of scatterers (3.97 X 10^^ 

N'/mL) was used for crude oil because of its ultra- 
absorptive characteristic. All lidar data are listed in 
Appendix B. 

Approximate values of the extinction coefficients 
for the different test oils were determined experimentally 
or from previous work [36] for laser wavelengths of 
.6328 urn and .4416 pm. Experimental values for the 
extinction coefficients were obtained by measuring the 
attenuation of laser light through a 1 mm wide glass sam- 
ple cell containing the different test oils and then 
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calculating the values from the exponential attenuation 
^ ot t 

formula e 
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CHAPTER VI 
RESULTS 

The previous chapters have investigated, analytically 
and experimentally, the problem of explaining the back- 
scatter from turbid water and oil on turbid water. Chapter 
III presented the development for backscatter from turbid 
water considering both single scatter and multiple scatter 
effects. The last section of Chapter III extended these 
results to include the effect on the backscatter of a 
layer of oil on the turbid water, considering the refrac- 
tive index change and attenuation due to the presence of 
oil. Chapters IV and V described the laboratory and field 
experiments which were conducted to: 1) determine the 

validity of the turbid water and oil on turbid water 
scattering models; 2) obtain a better understanding of the 
nature of the scattering phenomena; and 3) determine the 
applicability of a lidar polarimeter for remote measure- 
ments of subsurface water quality parameters and oil on 
water. This chapter describes the results of the investi- 
gation by analyzing the turbid water data, the natural 
waterway data, and the oil on water data, and comparing 
these results with the analytical models. 
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Turbid Water 

Chapter III presented the development o£ an analytical 
model for backscatter from smooth surfaced, turbid water 
by modeling the like polarized return as an incoherent sum 
of a single scatter term, (III-9), and a multiple scatter 
term, (111-13), where the single scatter term could be 
uniquely determined from the characteristics of the 
scattering medium. The analytical model for the cross 
polarized return consisted only of a multiple scatter term, 
represented by (III-14). This section first considers the 
validity of the single scatter term by examining the like- 
polarized returns from an aqueous suspension having low 
particle concentrations, where single scatter should be the 
dominant scattering mechanism. This section then considers 
the behavior of the entire analytical turbid water model 
(represented by (III-15) and (III-16) by comparing the 
relative magnitudes of the like-polarized volume reflection 
coefficients (y^^ and , and the cross polarized volume 

reflection coefficients (YyjJ) as a function of particle 
concentration. The results of this comparison will show 
that at low particle concentrations single scatter is the 
dominant scattering mechanism (y^^ is much greater than 
'^vh^ ’ and that at very high particle concentrations, 
multiple scatter is the dominant scattering mechanism 
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and Yy^ are much greater than Yyy) • 

Finally, this section discusses several characterise 
tics of measured backscatter that were observed from varia- 
tions of particle concentration, receiver field of view, 
transmitter polarization, and aqueous absorbent concentra- 
tion. 

Comparison With Analytical Model 

In examining the validity of the smooth surfaced, 
turbid water model developed in Chapter III, the subsurface 
single scatter term of (III-15) is first considered. Theo- 
retical values of single scatter return were computed for 
varying polystyrene latex particle concentrations and 
fields of view using (III-9). The system parameters used 
for these calculations are given in Table VI- 1 and Appendix 
B. The relative scattering function, F^(tt) , for the latex 
particles was calculated using an IBM 360 computer routine 
developed by Adams [37] . The relative volume extinction 
coefficient, 6^, was obtained experimentally by measuring 
the attenuation of laser light through very dilute concen- 
trations of latex particles (dilute concentrations were 
imperative for this parameter determination to insure 
measurement of single scatter attenuation only). The value 
of was not determined analytically because the imaginary 
part of the refractive index was not known and calculations 
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Table VI-1. Several Parameters Used for 

Comutations in Single Scatter 
Model 


Power Transmitted 


2 mw RMS 


Optical Efficiency 
Area of Receiver 


.5 

11.3 cm^ 


Table VI-2. Scattering and Extinction Coefficients for 
Polystyrene Latex and Teflon Particles 


Polystyrene Teflon 
Latex 

Relative Scattering 
Function (F (it)) 

(cm^/partitle) 1.033 x 10 1.163 x 10 


Relative Extinction 
Coefficient (6 ) 

(cm^/grams) 18219. 


81.9 
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by Kerker [33] show that the magnitude of the extinction 
cross section (or equivalently 3^.) is highly dependent upon 
the imaginary part of the refractive index. 

Theoretical values of single scatter return were also 
computed for an aqueous suspension having varying concen- 
trations of teflon particles. The value of F^(tt) was de- 
termined empirically by a data fit in the low range of 
particle concentration. The value of 3^. was determined by 
attenuation measurements, as discussed previously. The 
values of these coefficients for the latex particles and 
the teflon particles are given in Table VI-2. The relative 
scattering function and extinction coefficient for an 
average sized latex particle is several orders of magnitude 
greater than the coefficients for an average sized teflon 
particle. Figures V-1 (p. 54) and V-3 (p. 55) show that 
the mean size of the latex particles is several times 
larger than the mean size of the teflon particles. Since 
it is well known that larger particles scatter much more 
incident radiation than do smaller particles, the values 
for these coefficients appear reasonable. 

Employing these coefficients in (III-9), the theo- 
retical and experimental values of like polarized single 
scatter return were plotted as a function of latex parti- 
cle concentration (Figure VI-1) and teflon particle con- 
centration (Figure VI-2). Figures VI-1 and VI-2 indicate 



Particle Concentration (xlO particles/mL) 


Figure VI-1. Comparison of Calculated Single Scatt 
Return with Experimental Data as a 
Function of Field of View and Latex 
Particle Concentration 





Return Power (natiowatts) 



Figure VI-2. Comparison of Transmit Polarization 
and Calculated Single Scatter Return 
with Experimental Data as a Function 
of Teflon Particle Concentration 
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that at very low particle concentrations, the experimental 
backscatter is proportional to particle concentration, 
which is consistent with the expression for the single 
scatter volume reflection coefficient, given by 

CIII-11). Observe also from these figures that as the 
particle concentration becomes larger, the experimental 
backscatter tends to saturate briefly and follow the asym- 
tote of the calculated single scatter return. Finally, 
Figure VI -1 shows that at low particle concentrations , the 
experimental single scatter varies proportionally with the 
receiver field of view (or equivalently the mean depth L, 
given by (IV-1)), in direct agreement with the expression 
for the single scatter volume reflection coefficient, 
given by (III-ll) . 

The good agreement obtained in comparing the experi- 
mental values with calculated values (at particle concen- 
trations where single scatter appears to be the dominant 
mechanism) suggests that the single scatter term of 
(III-15) is valid. For large particle concentrations, 
however, the single scatter term does not adequately de- 
scribe the backscatter and, therefore, the entire analy- 
tical model (represented by (III-15) and (III-16)), which 
incorporates multiple scatter effects, must be used. The 
following paragraphs examine the behavior of the entire 
model. 
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The validity o£ the entire analytical model is best 
established by comparing the relative magnitudes o£ the 
single scatter volume reflection coefficient, with 

the multiple scatter volume reflection coefficients, 
and Referring to (III-15) and (III-16) , observe 

that the multiple scatter volume reflection coefficients, 
and Yy^> represent the fraction of incident radiation 
that is multiple scattered in the direction of the receiver. 
The effective single scattter volume reflection coefficient, 
Y^y, represents the fraction of incident radiation that is 
single scattered to the receiver. Examining the relative 
magnitudes of these volume reflection coefficients as a 
function of particle concentrations should provide a 
better understanding of the nature of the scattering 
mechanism and also establish the validity of the entire 
analytical model. 

The magnitudes of the volume reflection coefficients 
are plotted in Figure VI-3 as a function of latex particle 
concentration for a field of view of .64 cm. The curve 
representing Yyy was calculated from the relative scatter- 
ing coefficient, the mean depth L, and the particle con- 
centration, using (III-IO). The values of the like polar- 
ized multiple scatter coefficient y^y were determined by 
solving (III- 15) for y^^^ and using the measured values of 
like polarized backscatter for and the above mentioned 



Volume Reflection Coefficient (x 10 
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values for The values for Yy]J were calculated from 

(III-16) using the measured values of 

9 

For particle concentrations less that 10 particles/ 

s 

ml, the single scatter volume reflection, Yyy> is much 
greater than the cross polarized multiple scatter volume 
reflection coefficient Yy^* Therefore, in this particle 
concentration range, single scatter appears to be the 
dominant scattering mechanism (the values of Yyy at parti- 
cle concentrations where single scatter is dominant should 
not be used for comparison since these values represent the 
error in the single scatter model). However, observe that 
as the particle concentration increases beyond this region, 
the multiple scatter coefficients begin to become a more 
significant fraction of the backscattered energy. At a 

particle concentration of approximately 10^*^ particles/ml, 

s m 

Y^y appears insignificant compared to the values of Yyy 

or Yy^. that is, multiple scatter appears to be the domi- 
nant scattering mechanism. Observe that at even higher 
concentrations of particles, Yy^ approaches Yyy> which 
implies that the polarization of the backscattered wave is 
completely random and has no preferred orientation. 

Observing Figure VI-3 (p. 71) also shows that the 
cross polarized volume reflection coefficient increases 
steadily with particle concentration. The increase in' 

Yyh the single scatter region at concentrations less 
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than 10^ particles/ml are likely to be due to the minute 
nonspherical properties of the particles and the presence 
of a small amount of multiple scattered radiation. At 
concentrations much higher than this, the increase in 

with particle concentration is due to increased multi- 
ple scatter. 

The previous examination of the relative magnitudes 
of the volume reflection coefficients shows that each 
volume reflection coefficient (Yyyi ■Yvh^ repre- 

sents a significant fraction of the backscattered energy 
during certain portions of the progression from low 
particle concentrations to high particle concentrations. 
Therefore, the entire analytical model for describing the 
backscatter from turbid water, represented by (111-15) 
and (Ill-lb), appears valid. Previous analytical back- 
scatter models such as those developed by Rouse [30] , 
Leader [29], and Wilhelmi [18] have ignored the presence 
of single scatter from the subsurface. It is apparent 
from the results presented in Figures VI-l(p.67), VI-2 
(p. 68), and VI-3 (p. 71) that these effects cannot be 
ignored, except for a very dense medium, where multiple 
scatter is the dominant scattering mechanism. 
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Discussion 

Having determined the validity of the analytical model, 
a discussion of the characteristics of the data described 
in Chapter V will be given to better understand the nature 
of the scattering phenomena. The first effect examined is 
the relationship between the presence of significant mul- 
tiple scatter and the particle concentration. The calcu- 
lated single scatter return plotted in Figure VI-1 (p* 67) 
appears to be consistent with experimental data when the 
value of the argument of the exponential in (III-9), 

(2N 8^L) , is less than approximately 1.2, and appears to 
vary proportionally with concentration for arguments less 
than approximately 0.2. These attenuation values are in 
fair agreement with the guidelines set forth by van de 
Hulst [32] in giving the critical attenuation values for 
determining the presence of multiple scatter. Van de Hulst 
claims that whenever the value of the argument in (III-9) 
is greater than 0.6, multiple scattering effects should be 
considered, and that single scatter should be the dominant 
mechanism for values of the argument less than 0.2. 

While the attenuation factors mentioned above are 
reasonable indicators for determining when multiple 
scattering becomes dominant, possibly a stronger indication 
of the occurrence of this phenomena is the depolarization 



ratio (cross polarized return divided by like polarized 
return). Figures VI-4 and VI-5 show that the depolariza- 
tion ratio is small and relatively constant when single 
scatter is the dominant mechanism. Observe that it is 
when the single scatter term saturates that the depolar- 
ization ratio begins to increase, that is, when multiple 
scatter becomes the dominant scattering mechanism. 

An interesting comparison can be made concerning the 
different depolarization ratios that are observed from the 
backscatter from low concentrations of latex particles 
and teflon particles. Observe that the depolarization 
ratio value for low concentrations of teflon particles 
(Figure VI-5) remains essentially constant at 0.1, while 
the depolarization ratio for the latex particles (Figure 
VI-4) remains essentially constant at a much smaller value 
of 0.03. The difference in the depolarization ratio value 
can be explained by observing Figures V-2 (p. 54) and 
V-4 (p. 55). Observe that the sphericity of the latex 
particles is much greater than that of the teflon parti- 
cles. The latex particles, therefore, depolarize the 
incident radiation much less that the teflon particles 
(it is generally accepted that nonspherical scatterers 
depolarize incident radiation). Since this phenomena 
should introduce depolarization proportional to particle 
concentration, the depolarization ratio at low concentra- 



Depolarization Ratio 
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Figure VI -5. Graph of Depolarization Ratio as a 
Function of Teflon Particle Con- 
centration ■ 



78 


tions of scatterers (where single scatter is proportional 
to concentration) should remain constant. 

Another effect studied to better understand the 
nature of the scattering mechanism was the effect on the 
backscatter of varying the field of view of the receiver. 
From Figures VI-1 (p. 67), VI-4 (p. 76), and VI-6, it is 
observed that not only do the like and cross polarized 
backscatter returns increase with increases in field of 
view, but the depolarization ratio increases as well. This 
is significant in that it proves that the depolarization 
is a multiple volume scatter process. This is apparent 
from the fact that if the depolarization was due to the 
nonspherical properties of the particles, the like- 
polarized and cross-polarized returns would increase 
proportionally with an increase in volume. However, since 
the cross-polarized returns increase more than the like 
returns (due to an increase in volume) , the depolarization 
process must be due to scattering outside of the incident 
beam, that is, a volume scattering process. 

Another effect examined in the turbid water experi- 
ments was the effect on the backscatter of adding an absor- 
bent to the scattering medium. From Figure VI-7, it is 
observed that the depolarization ratio varies inversely 
with absorbent concentration. This observation supports 
the presence of a multiple scatter mechanism since the 
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Graph o£ Depolarized Return as 
a Function of Receiver Field of 
View and Latex Particle Concentration 


Figure VI -6. 
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Figure VI-7. Graph of Depolarization Ratio as 

a Function of Absorbent Concentration 
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total path length of a multiple scattered wave is longer 
than that of a single scattered wave. The exponential 
attenuation Ce’^^, where x is the path length) due to the 
absorbent, decreases the intensity of the multiple scat- 
tered wave much more than the single scattered wave. This 
effect could hinder measurements of the particle concen- 
tration, since the behavior of the absorbent effects the 
depolarization ratio in the same manner as would the effect 
of decreasing the particle concentration. Note, however, 
that large dye concentrations were required to vary the. 
depolarization ratio significantly and that under normal 
conditions in the natural environment it is not expected 
that the absorbent concentration would vary as dramatically. 

The last effect studied in the turbid water experiment 
was the effect on the backscatter of changes in the trans- 
mit polarization. As can be observed from Figure VI- 2 
(p. 68) , the like polarized returns were fairly consistent 
and independent of transmit polarization. The cross- 
polarized returns for horizontal transmit polarization 
appear to be a small amount larger than the cross-polarized 
returns for vertical transmit polarization (the cross 
returns at very low particle concentration of particles are 
beyond the noise limitation of the system and should not 
be used for comparison). This small inconsistency of 
cross returns does not support the analytical model, that 
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is, the analytical model supports the law of reciprocity. 
Other investigators such as Rouse [30] also contend that 
recriprocity must hold for an isotropic media. Therefore, 
the small differences of and must be due to the 
presence of nonspherical scatterers. 

Natural Waterways 

Representative results of the first field measurement 
trip (data are listed in Appendix B and coded by 71473 and 
71573 under DATE) from atop a bridge on the Brazos River 
near Waco, Texas are plotted in Figures VI-8 to VI-11, 
Relatively few data observations were made during this 
trip because of several experimental difficulties. These 
data, however, did provide initial indications that the 
like polarized backscatter, cross polarized backscatter, 
and the depolarization ratio vary proportionally with 
turbidity and suspended solids and inversely with trans- 
mittance. The twenty minute time history of P^^ (Figure 
VI- 10) shows clearly that an increase in turbidity results 
in an increase in Pyy* Figure VI-10 also implies that for 
that particular location, the ordinate of the time history 
could have essentially been turbidity and provided a 
continuous, real-time measurement of that parameter. The 
few horizontal transmit polarization measurements showed 
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Figure VI- 8. Graph o£ Lidar Returns as a Function o£ 
Suspended Solids Observed £rom a Bridge 
on the Brazos River 
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Figure VI -9. Graph of Lidar Depolarization Ratios as a 
Function of Suspended Solids Observed from 
a Bridge on the Brazos River 
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Figure VI -11. Graph of Lidar Depolarization Ratio 

as a Function of Turbidity Observed 
from a Bridge on the Brazos River 
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Figure VI -12. Graph of Like Polarized Returns as a 

Function of Transmittance Observed along 
Houston Ship Channel from a Boat 
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insignificant differences in backscatter relative to ver- 
tical transmit polarization measurements. 

A representative plot of the results from the field 
measurement trip on the boat Excellence along the Houston 
Ship Channel (data are listed in Appendix B and coded by 
92573 under DATE) is given in Figure VI-12 (p. 85). Note 
that the return power supports the results of the previous 
field measurement trip by varying inversely with trans- 
mittance. A significant measurement error was encountered 
in these measurements which partially explains the scat- 
tered data points. The lidar slant range aboard the 
Excellence was approximately 2.84 m, with a receiver field 
of view diameter at the surface of 1.26 cm and a distance 
between the receiver and transmitter of 12.4 cm. According 
to (IV-1), the values correspond to a mean depth along the 
beam of approximately 14.4 cm. Since the peak value of 
water wave heights was approximately 10 cm, the lidar 
measurement volume was changing significantly with wave 
height variations. Time averaging with the lock-in ampli- 
fier smoothed these variations in returns, but did not 
provide accurate and consistent measurements of back- 
scatter. As a result of this problem, a system constraint 
was developed for future field measurements to improve 
performance in a natural environment with water waves. 

This constraint requires that the mean depth, L, multiplied 
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by cos 6 (since L is the depth along the incident beam) 
be much greater than the peak wave height^ A, or 


L cos(e) » ^ 


which when substituting (IV- 2) for L' gives 


R^co s(e) 
vTf 


» 2 A 


(VI-l) 


This equation represents a parametric decision criterion 
for lidar operation under field conditions. To minimize 
wave effects using this equation requires careful consid- 
eration because several trade-offs are involved. Increas- 
ing the receiver field of view diameter (d^j^R/F) does tend 
to minimize wave effects, but only at the expense of 
increased background noise (solar radiation) . An optimal 
procedure then for the minimization of wave effects is 
first to minimize the distance, w, between the transmitter 
and receiver. Second, for a given altitude, A, the slant 
range, R, is given by 


R 


A /cos( e) 
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which when substituted in (VI- 1) gives 


A^d 


ph 


w cos(e)F 


» 2 A 


This equation shows that wave effects can also be minimized 
by operating at grazing angles. However, observations of 
the Fresnel transmission coefficient show that these coef- 
ficients vary dramatically with changes in incidence angles 
near grazing. Wave variations, that is, local angle varia- 
tions, would therefore create significant backscatter 
’’noise” at incidence angles greater than 50". A general 
procedure for minimizing wave effects considering all pa- 
rameters is to: 1) minimize the distance between the 

transmitter and receiver; 2) operate at incidence angles 
less than 50° but away from nadir (to avoid specular re- 
flections) ; and 3) maximize the slant range without forcing 
field of view at the surface to be much larger than the 
diameter of the incident beam. Observe that a receiver 
field of view matched to the transmitted beam is not opti- 
mum, since wave variations can significantly effect the 
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measurement volume. Therefore, at ranges less than about 
20 m, the diameter of the receiver field of view at the 
surface should be at least twice the diame-ter of the inci- 
dent beam. 

Employing these criteria on the last field measurement 
trip aboard the Excellence (data are listed .ia Appendix B 
and coded by 120173 under DATE) provided excellent results. 
Representative results are plotted in Figures VI-13 to 
VI-15. Note that the data consistency is better than the 
data from previous field measurement trips and that the 
depolarization varies significantly with changes in tur- 
bidity, suspended solids, and transmittance. Also the 
magnitude of the depolarization ratios are within the same 
range as those measured from simulated turbid water in the 
laboratory. Therefore, the laboratory simulations of 
turbid water employed earlier to test the validity of the 
analytical model are valid simulations of natural water. 

The magnitudes of the depolarization ratios are also within 
the range of laboratory depolarization ratios which varied 
directly with particle concentrations. This is significant 
because this confirms that the magnitudes of particle con- 
centrations that are found in natural waterways are large 
enough to cause changes in depolarization ratios for dif- 
ferent concentrations. These values also infer that single 
scatter is also present from the subsurface, and that the 
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Figure VI -14. Graph of Lidar Depolarization Ratios as 

a Function of Turbidity Observed Along 
Brazos River From a Boat 
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Figure VI-15. Graph of Lidar Depolarization 

Ratios as a Function of Trans- 
mittance Observed Along Brazos 
River From a Boat 
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model found valid in the previous sections appears valid 
for describing the backscatter from natural waterways. 

An experimental error common to all field measure- 
ment data was the requirement of a 30 second time interval 
between consecutive recordings of and due to the 

single channel capacity of the lock-in amplifier. This 
time lag between recordings coupled with collecting 
water samples within these intervals can add significant 
error (approximately 5%) because of the dynamic changes 
that occur in the subsurface over a short period of time, 
as shown in Figure VI-10 (p. 84). 

Oil on Water 

Comparison Vifith Analytical Model 

The analytical model for oil on water represented by 
(III-20) and (III-21) incorporated the effects on the back- 
scatter of the change in refractive index and attenuation 
due to the presence of a laterally inhomogeneous layer of 
oil on water. It was assumed in this model that the at- 
tenuation due to oil was independent of polarization, that 
is, the oil attenuated both like and cross polarized re- 
turns in the same fashion. Employing these results, the 
theoretical values of backscatter were calculated for 
varying thicknesses of different test oils on water. 
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Representative values for the indices of refraction, trans- 
mission coefficients, and homogeneous extinction coeffi- 
cients employed in (111-20) and (III-21) were determined 
empirically or from published work [36] , and are given 
in Tables VI- 3 and IV-4 (the significance of these values 
is discussed in the next section.) For these laboratory 
studies, the oil layers were laterally homogeneous, which 
allowed use of homogeneous indices of refractions and 
extinction coefficients (p=l in (III-19) and (III-22)). 

The value for and in (III-20) and (III-21) are 
the measured values of backscatter in the absence of oil. 

Graphical comparisons of calculated values and ex- 
perimental values of backscatter for several test oils 
(values are for an excitation wavelength of X=,6328 ym) 
are given in Figures VI-16 to VI-19, as a function of oil 
thickness. Observe from Figures VI-16 and VI-17 that the 
attenuation due to crude oil is quite significant, re- 
sulting in a decrease in like or cross polarized back- 
scatter greater than 3 orders of magnitude for an oil 
thickness of .3 mm. The attenuation due to no. 2 fuel 
oil is less than that of crude oil, but still results in a 
decrease of backscatter, greater than 1 order of magnitude 
for a 1 mm oil thickness. 
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Table VI-3. Values of Refractive Indices and Transmission 
Coefficients Used in Calculations of Back- 
scatter from Oil on Water 


Test 

Oils 

T T V 

T\ 1-1 ly fc i « -| 

0 13^ 31 

V _ V 
12* '21 

.p V „ V 

''3 2 ’ ^2 3 

T. h 

'32 

rp h rp h 

21 31 

Crude 

Oil 

1.48 

.976 

.997 

.996 

.945 

No. 2 
Fuel Oil 

1.44 

.980 

.999 

.998 

.953 

SAE 30 

1.43 

.981 

.999 

.998 

.954 

Gasoline 

1.37 

.985 

.999 

.999 

.963 

Kerosene 

1.37 

.985 

.999 

.999 

.963 

Water 

1.33 .988 

— 



— 

.970 


Table VI-4. Homogeneous Extinction Coefficients for the 
Test Oils at Two Different Wavelengths 



a 


Ratio 

Test 

(nepers /mm; 

(nepers/mm; 

(a /a. 

Oils 

X=.6328 ym) 

X-.4416 ym) 

1 D 

Crude 

Oil 

12.5 

9.12 

1.37 

No. 2 

Fuel Oil 

1.5 

7.10 

.207 

SAE 30 

.05 

.35 

. 143 

Gasoline 

.006 

.07 

.08 

Kerosene 

.015 

no measurable 
attenuation 

— 
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Figure VI-16. Comparison of Analytical Model with 

Experimental Data for Like Polarized 
Backscatter from Crude Oil and No. 2 
Fuel Oil on Turbid Water 
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Figure VI-17. Comparison of Analytical Model with 

Experimental Data for Cross Polarized 
Backs catter from Crude Oil and No. 2 
Fuel Oil on Turbid Water 
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Figure VI- 18. Comparison of Analytical Model with 

Experimental Data for Like Polarized 
Backs catter from Kerosene on Water 
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Figure VI- 19. Comparison of Analytical Model with Experi- 
mental Data for Like Polarized Backscatter 
from SAE 30 and Gasoline on Turbid Water 
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The attenuation due to kerosene is minimal (Figures 
VI-18 (p. 98) and VI-19 (p. 98)), even at oil thicknesses 
as large as 6 mm. The attenuation due to gasoline (Figure 
VI-19 (p. 98)) is also quite small for the red excitation 
wavelength used in these measurements, decreasing the 
backscatter only a few percent for large oil thicknesses. 
The attenuation due to the presence of SAE 30 refined motor 
oil is also relatively small, decreasing the backscatter 
approximately 15 percent for an oil thickness of 1 mm. 

Good agreement was obtained between calculated values and 
experimental values of backscatter for all test oils 
(Figures VI-16 to VI-19 (pp. 96-98)). An examination of 
the depolarization ratio will now be given to determine 
its behavior as a function of oil thickness. 

The depolarization ratio remains essentially constant 
for different thicknesses of oil on water (Figures VI- 20 
to VI-22). This is in direct agreement with the calculated 
values for the depolarization ratio, since the model as- 
sumed equal attenuation for like and cross polarized re- 
turns. Large values of depolarization ratios were observed 
because of the large particle concentrations that had to 
be used (for large thicknesses of crude oil) to provide 
magnitudes of backscatter above the noise limitations of 
the system. Therefore, any volume scatter from the oil 
itself may have been masked by the large contribution of 
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Figure VI -20. Comparison of Analytical Model with 

Experimental Data for the Depolar- 
ization Ratios from Crude Oil 
on Turbid Water 
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Figure VI-21. Comparison of Analytical Model with 
Experimental Data for the Depolar- 
ization Ratios from No. 2 Fuel Oil. 
and SAE 30 on Turbid Water 
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Figure VI- 22. Comparison of Analytical Model with 
Experimental Data for the Depolar- 
ization Ratios from Gasoline and 
Kerosene on Turbid Water 
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volume scatter from the subsurface. It is felt, however, 
that the volume scatter contribution from any of the test 
oils would be much smaller than that observed from the 
subsurface . 

The excellent agreement between the calculated values 
and experimental values of backscatter and depolarization 
ratios suggests that the analytical model for oil on water 
represented by (III-20) and (III-21) is valid. 

The estimated experimental error in the data acqui- 
sition of the oil on water experiments is 5 percent. The 
most significant error which was encountered in the pro- 
cedure was the experimental measurements of oil thickness. 
Oil thickness was measured by adding a known volume of oil 
to a known, contained surface area, with the quotient of 
the volume and surface area giving the oil thickness. The 
known volume of oil added was subject to several experi- 
ment constraints such as the adhesive characteristics of 
the oil to the measurement vessle. The oil on water also 
tended to move toward the edge of the container which re- 
duced the oil thickness at the center of the container, 
where the laser beam was incident. 
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Discussion 

The detection, identification, and measurement feasi- 
bility of oil on water can be established by observing the 
following characteristics of the data presented in Figures 
VI-4 Cp. 76) and VI-5 (p. 77), VI-16 to VI-19 (pp. 96-98), 
VI-20 to VI-22 (pp. 100-101) , and in Table VI-4 (p. 95); 

1) the variation of subsurface particle concentrations 
significantly alters the depolarization ratio; 2) the 
presence of oil does not significantly alter the depolar- 
ization ratio; and 3) the attenuation due to oil is a 
function of oil type, oil thickness, and excitation wave- 
length. Using these results, a general procedure can be 
followed for detection, type identification, and thickness 
measurement of oil on water. The detection of oil on water 
should require only a single wavelength lidar polarimeter 
and an observance of and the depolarization ratio 

(Pyh-^^vv^ ‘ significant changes in P^^ or without 

a significant change in depolarization ratio should indi- 
cate the presence of oil. 

Oil type identification can be determined using a 
lidar polarimeter by incorporating additional spectral 
information in the backscatter. Employing a dual wave- 
length lidar polarimeter system such as that being con- 
structed by the Remote Sensing Center at Texas A§M 
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University should provide adequate spectral information to 
perform oil type identification. Referring to the analy- 
tical model for oil on water given by (III-20), observe 
that the extinction coefficient in (III-20) has been shown 
from Table VI -4 (p. 9 5) to be a function of oil type and 
excitation wavelength. Therefore, the like and cross 
polarized returns measured by a dual wavelength system can 
be represented by the following equations: 


^vv 

= T'exp(-2.,t) P^; 

(VI-2) 

r 

Pvh 

= T"exp(-2a, t) 

(VI-3) 

•vv 

= T'exp(-2«,t) Py'; 

(VI-4) 

pb 

Vh 

= T"exp(-2.,t)P^;; 

(VI-5) 


where designates the red returns in the absence of oil; 
P^j designates the blue returns in the absence of oil; and 

and T’ and T" represent the transmission coefficients in 

\ 

(III-20) and (II-I-21), Oil detection is first established 
by following the same procedure mentioned above for the 
single wavelength system. Once oil detection is estab- 
lished, the ratio of the natural logarithm of P^^ and P^^ 
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should be observed from which the ratio can be cal- 

culated as 



Referring to Table VI-4 (p. 94), observe that this ratio is 
single valued for the oil types, that is, each extinction 
coefficient ratio corresponds to only one oil type. There- • 
fore, employing this ratio as an indicator of oil type, 
allows identification of oil type through dual wavelength 
backscatter measurements. 

Oil thickness determination follows easily after oil 
type has been established by employing a table look-up 
procedure to determine the value of the extinction coef- 
ficient for the specific oil type. Employing one of these 
extinction coefficients, for instance a^, in (VI-2) allows 
solution for the oil thickness t (ignoring transmission 
coefficient effects) . 

A possible ambiguity involved in the above procedure 
is that the extinction coefficient may not correspond to 
only one oil type. However, for the test oils used in this 
study, this ambiguity was not present. The detection 
scheme ignores the transmission coefficient effects due to 
a layer of oil on water. Observing the magnitudes of the 
transmission coefficients listed in Table VI-3 (p. 94), 
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however, shows that the product of the transmission coef- 
ficients used in (III-20) and (III-21) is greater than .95. 
Therefore, if these effects are ignored, a maximum error of 
5 percent in calculating will be encountered. Smaller 
errors can be obtained by including a table look-up pro- 
cedure for the index of refraction, once oil type is de- 
termined. This error does not significantly effect oil 
type identification (since ratios are used) , but does af- 
fect oil thickness determination. 

Table VI-4 (p. 94) and Figure VI-18 (p. 97) show 
that the attenuation of laser light due to kerosene is 
minimal. Therefore, kerosene on water does not appear to 
be detectable by the wavelengths observed in this study. 
However, as reported by Horvath et a^. [36], kerosene 
has a strong, narrow absorption band at a wavelength of 
approximately .3 pm and could therefore be detected by a 
lidar polarimeter extending its wavelength to the ultra- 
violet region of the spectrum. 
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CHAPTER VII 

CONCLUSIONS AND RECOMMENDATIONS 

This investigation has examined, experimentally and 
analytically, the applicability o£ a dual polarization 
laser backscatter system for remote identification and 
measurement of subsurface water turbidity and oil on water. 
The following paragraphs describe the conclusions drawn 
from this study and the recommendations for future work. 

Conclusions 

The analytical model for describing the backscatter 
from smooth surfaced, turbid water developed in this in- 
vestigation included single scatter and multiple scatter 
effects. A comparison of the experimental data with the 
analytical model showed that the analytical model was valid 
and that single scatter effects cannot be ignored, except 
for scattering from a very dense medium. The experimental 
laboratory data for backscatter from turbid water also 
showed that the depolarization ratio: 1) remained constant 

as a function of particle concentration at low particle 
concentrations; 2) varied directly with particle concen- 
tration at high particle concentrations; 3) varied directly 
with the receiver field of view; and 4) varied inversely 
with absorbent concentration. Experimental laboratory 
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data also showed that changes in transmit polarization had 
an insignificant effect on the backscatter . 

Field measurement data from several natural waterways 
showed that the like and cross polarized returns and the 
depolarization ratio varied directly with the water 
quality parameters turbidity and suspended solids, and 
inversely with transmittance. These results support the 
use of a lidar polarimeter for remote measurements of tur- 
bidity, suspended solids, and transmittance. 

The analytical model developed in this investigation 
for describing the backscatter from oil on water included 
the effect on the backscatter of the change in refractive 
index and attenuation due to the presence of oil on water. 
Experimental laboratory data with laterally homogeneous oil 
layers verified the validity of the analytical model and 
established the following characteristics of the data: 

1) the presence of oil on turbid water does not signifi- 
cantly alter the depolarization ratio; and 2) the attenua- 
tion due to oil is a function of oil type, oil thickness, 
and excitation wavelength. Considering these effects on 
the backscatter, a detection scheme was presented whereby 
the presence of oil on water can be established using a 
single wavelength lidar polarimeter. A detection, type 
identification, and thickness measurement scheme was also 
presented using a dual wavelength lidar polarimeter. These 
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results strongly support the use of a lidar polarimeter for 
oil pollution detection and measurement. 

Recommendations 

Based on the findings of this study, several recom- 
mendations are made concerning the development of a more 
exact model of the backscatter from turbid water and the 
enhancement of the applicability of a lidar polarimeter 
for remote measurements of water pollution. The like 
polarized single scatter model presented in this investi- 
gation included an exact treatment of the single scattering 
mechanism, that is, the single scatter model was given in 
terms of the physical properties of the medium, and, there- 
fore, required empirical results to determine its behavior. 
Future investigations should include an analytical treat- 
ment of the multiple scatter volume reflection coefficients 
(cross and like polarized) to derive expressions for these 
coefficients as a function of one another or as a function 
of the medium properties. 

To more firmly establish the applicability of a lidar 
polarimeter for remote measurements of water quality param- 
eters, a more exhaustive field measurements program should 
be conducted to determine exact relationships between the 
lidar returns and the physical water quality parameters. 
Spectral characteristics of backscatter from natural 



110 


waterways should also be investigated for possible identi- 
fication of harmful, chemical subsurface pollutants. 

Future investigations for oil pollution detection and 
identification should consider more thoroughly the spectral 
characteristics of many oils. This study is necessary so 
that the assumptions employed in the identification scheme 
(the ratio of two extinction coefficients for two different 
wavelengths is single valued) can be more thoroughly 
tested. Also, future experimental investigations should 
determine if the inhomogeneities of the oil layer are cor- 
rectly accounted for in the analytical model, so that these 
effects may be incorporated into the detection scheme. 
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APPENDIX A 

DESCRIPTION AND DIAGRAM OF LIDAR 
POLARIMETER CONTROL UNIT 

The most significant improvement of the second 
generation field portable lidar polarimeter was the 
packaging of a central electronic control unit (see Figure 
A-1) which allows convenient operation, field portability, 
and faster experimental set-up time. This appendix des- 
cribes a few of the considerations that were involved in 
the design of the 'electronicnunit and gives detailed 
schematics and wiring diagrams of the electronic circuitry. 

Problem Description ' ’ 

There are several electronic, components contained 
within the lidar polarimeter: / laser, >laser pbw supply, 
two photodiode/operational amplifier (op-amp) combination 
detectors, photodiode/op- amp power supply, light chopper, 
and lock-in amplifier. The signals transferred between 
the above components include: chopping reference frequency 

light level signals from each of the two detectors, DC 
offset bias voltages for each op-amp, and driving power 
inputs for all components. 
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Figure A- 1 . Lidar Polarimeter Electronic 

Control Unit 


Design Considerations 


The magnitude of the light level voltages trans- 
ferred from the photodiode op-amp detectors can be on the 
order of tens of microvolts. Noise design considerations 
such as proper shielding of signal wires and operating com- 
mon mode were imperative. 

Another consideration involved in the design of the 
electronic unit was the requirement of a twenty- five 
minute warm-up period for the laser power output to reach 
its steady state peak value. Use of the lidar polarimeter 
during this warm-up period would yield inaccurate measure- 
ments. A special timing circuit was therefore implemented 
to notify the lidar operator when this twenty- five minute 
warm-up period expired. Digital gating circuitry for all 
components was also used to notify the operator when all 
system components are activated and ready for operation. 

A final consideration involved in the design of the 
electronic unit involved adapting the lidar system to the 
power supplied by a gasoline-powered 110 VRMS, 60-cycle 
generator needed for field operation. Field tests with 
the lidar system powered by a gasoline generator showed 
that the frequency stability of the 1500 watt generator 
was inadequate. The light chopper had a synchronous motor 
and the chopping frequency would not remain sufficiently 
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stable to allow synchronous detection with the lock-in 
amplifier. To alleviate this problem a special push-pull 
power amplifier circuit was designed and constructed to 
amplify the voltage supplied by a stable 60-cycle generator 
[Model 120-Exact Generator] to the necessary voltage and 
current required by the light chopper. Field tests with 
this additional circuit resulted in successful operation 
of the lidar system when powered by a gasoline-powered 
generator . 


Electronic Unit Description 

The interconnection of the electronic control unit 
with the transceiver of the lidar polarimeter is diagrammed 
in Figure A- 2. Figure A- 3 gives the connector and wire 
list, listing all labeled connectors and wires with their 
manufacturer and type. Figure A-4 shows the detailed 
wiring diagram for the connections indicated in Figure A- 2. 

A schematic drawing of the 110 VRMS power wiring 
connections of the electronic control unit is given in 
Figure A- 5. The wiring of the indicator lights and other 
DC powered components is shown in Figure A- 6 with an 
assembly drawing for Figures A- 5 and A- 6 given in Figure 
A- 7. The timing, gating, and lamp driver circuits are 
diagrammed in Figure A- 8. 




Figure A-2. Lidar Polarimeter Electronic Connections, 




1J1,2J1,5J1 - Winchester 7 pin» M7S 
1P1,2P1,5P1 - Winchester 7 pin, M7P 
3J1 - Kings 79-46 

3Pl,5Pll-13,6Pl-3 - Amphenol UG88/U 

3J2 - Cinch Jones P303-CCT-K 

3P2 - Cinch-Jones P303-CCT-L 

4P1 - Crimp Terminal 

5J2 - Winchester 10 pin, M710S 

5P2 - Winchester 10 pin, M710P 

5J3 - Winchester 4 pin, M4S 

5P3 - Winchester 4 pin, M4P 

5J4 - Kings KV-79-18 

5P4 - Kings KV-59-27 

5J5-5J9 - Amphenol 160-4-N 

5J6,5J7,5J8 - Terminal Plugs 

5J10 - Amphenol 160-5-N 

5J11-13 - Amphenol UG657/U 

6J1-2 - Kings 79-46 

6J3 - Kings 79-35 

W(1),W(2) - Belden Audio, 3 shielded pairs. Type 8767 
W(3) - Belden Audio, 1 shielded pair, Type 8434 
W(4) - Belden H.V. COAX, RG/58U 
W(5),W(10) - llOV AC three conductor power cord 
W(6) - Single Conductor Insulated Wire (Gauge-22) 

W(7) - Belden COAX, RG158/U 


Figure A-3. Connector and Wire Parts List. 
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5J1 -A->5P1-A^-1P1 -A->lJl-A->-Pin 6, Op-Amp, signal 
5Jl-B->5Pl-B^lPl-B->lJl-B^Pin 7, Op-Amp, +15 VDC 
5J1-C-+5P1-C-+1P1 -C->-lJl-C-+Pin 5, Op-Amp, DC Bias 
5Jl-D-+5Pl-D-+lPl-D-+lJl-D-^Pin 1, Op-Amp, DC Bias 
5Jl-E^5Pl-E^lPl-E-+lJl-E-+Pin 3,4, ground 
5J2-A-»-5P2-A-+2Pl-A-+2Jl-A-+Pin 6, Op-Amp, signal 
5J2-B-v5P2,-B-^2Pl-B^2Jl-B+-Pin 7, Op-Amp, +15 VDC 
5J2-C-»-59 2-C-+2Pl-C-+2Jl-C-+Pin 5, Op-Amp, DC Bias 
5J2-D+-5P2-D'+2Pl-D->-2Jl-D-+Pin 1, Op-Amp, DC Bias 
5J2-E-+5P2-E-+2Pl-E-»-2Jl-E-+Pin 3,4, Op-Amp, ground 
5J3-^5P3-+3P2-+3J2'+cliopper reference 

5J4, switched high voltage+‘5P4 , laser high voltage 
SJ5, switched 110 VAC-+3P2 3J2, chopper 110 VAC 
5J6, Auxiliary + 15 VDC-+ 

5 J7 , Auxiliary Chassis ground->4Pl , optical base 
5J8, Auxiliary - 15 VDC+- 

5J9, Auxiliary 110 VAC-+5P9, lockin x 110 VAC 
5J10, 110 VAC input-^SJlO 
5J1-A-+5J11-+5P11^6P1^6J1, signal 
5J2-A->5J12^5P12-+6P2-»^6J2, signal 
5J3-A+-5J13-+5P13-+6P3-+6J3 , chopper reference 


Figure A- 4. Interconnection Diagram for 
Figure A- 2. 









5J1 5J2 5J3 5J4 5J5 5J6 5J7 5J8 5J9 5J10 



Control 
Unit Rear 
Panel 


Figure A- 6. DC Wiring Diagram, 
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PSl 

PS2 

PS3 

LI 

51 

52 

53 

54 

T1 

Cl, 

C2 


Spectra Physics Model 124 Laser Power Supply 

Wanlass OEM Power Supply, 5V @ 3A, Type DPSl-5 

Analog Devices Power Supply, ±15V, Model 950 

System Status Indicator Light, Dialco Type 533-0901, 
with lense Dialco Type 303-3472 

Main Power Switch, Dialco Type 513-1504-001 with 
lense Dialco Type 303-3474 

Laser Power Switch, Dialco Type 513-1509-001 with 
lense Dialco Type 303-3471 

Photo-diode OpAmp Power Supply Switch, Dialco Type 
513-1509-001 with lense Dialco Type 303-3473 

Chopper Power Switch, Dialco Type 513-1509-001 with 
lense Dialco Type 303-3475 

Cinch Terminal Strip - Type 6-141 

Spectrol lOK Precision Potentiometers - Model 1532 


Figure A- 7. Assembly Drawing for Figures A- 5 and A- 6. 



Ql = Q2=Q3=Q4=TI202 

Q 5 = 2N2369A 
A = SN7404 
B = SN7420 
C = NESS 5 



Timer, and Lamp Driver Circuits 
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The push-pull power amplifier circuit designed and 
constructed to provide a stable power frequency for the 
light chopper is diagrammed in Figure A- 9. Supplying a 
60-cycle 1.5 volt peak-to-peak voltage at the input of this 
circuit (along with a 35 V @ 7 amp operating voltage) 
results in a 110 VRMS output capable of supplying 14 watts 
of power. 



35 (3 7A 35V @ 7A 



Figure A-9. Push-Pull Amplifier Circuit. 


= Q 2 = 2N3439 
Q 3 » 2N404 
= 2N3442 
= 2N1046 
Di = D 2 = IN3605 
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APPENDIX B 
DATA 

This appendix presents the data collected for this 
investigation. The data is segmented into three types of 
measurements: 1) laboratory turbid water measurements 

(Table B-1); 2) laboratory oil on water measurements 
(Table B-2); and 3) field measurements from natural water- 
ways (Table B-3). Several lidar system parameters and 
measurement parameters that were used for this investi- 
gation are listed under the column CODE. The interpreta- 
tion for this code is discussed below. 

abed- The digit a represents the type of particles 
used for the turbid medium. Values for a are: 

1 - Dow Plastic Pigment 722 

2 - E.I. Dupont Resin Dispersion, TFE 30 

3 - Natural Water (field measurement) 

The digit b represents the slant range. Values for b are: 

1 - 9.5 m 

2 - 2.84 m 

3 - 1.44 m 

The digit c represents the diameter of the receiver field 
of view at the surface of the water. Values for c are: 
1-4.2 cm 
2 - 7.04 cm 
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3 - 21.1 cm 

4 - 70.4 cm 

5 - 1 . 26 cm 

6 - .32 cm 

7 - .64 cm 

8 ' 1.07 cm 

The digit d represents the distance between the 
and the transmitter at the transceiver section, 
for d are: 

1 - 12.4 cm 

2 - 7 . 3 cm 


receiver 

Values 



PARTICLE 

SAMPLE CONCENTRATION 


NUMBER 

DATE 

CODE 

(NUMflE«/ML ) 

PVVCWATTSJ 

PVHt WATTS! 

2 

11474 

1362 

0.175E 

07 

O.llTE-lO 

0.117E--12 

3 

11474 

1362 

0.354E 

07 

0.235E-10 

0.352E-12 

4 

11474 

1362 

0.530E 

07 

0.352E-10 

0.470F-12 

5 

H474 

1362 

0.123E 

08 

0.764E-10 

0.106E-11 

6 

11474 

1362 

0.192E 

OB 

O.U6E-09 

0.176E-11 

7 

11474 

1362 

0.367E 

08 

0.210E-09 

0.329E-11 

8 

11474 

1362 

0.715E 

08 

0.351E-09 

0.5876-11 

9 

11474 

1362 

0. 141E 

09 

0.603E-09 

0.1156-10 

10 

11474 

1362 

0.296E 

09 

O.033F-O9 

0. 186E-10 

11 

11474 

1362 

0-567E 

09 

o*i04E-oa 

0.322E-10 

12 

11474 

1362 

0- 142E 

10 

0*15lF-'08 

0.9O2E-10 

13 

11474 

1362 

0^272E 

10 

0 * 2341-08 

0.270E-09 

14 

11474 

1362 

0.735E 

10 

0.433b-0a 

0.1136-08 

15 

11474 

1362 

0.U7E 

ll 

a*62lC-08 

0.2356-08 

16 

11474 

1362 

0.160E 

11 

0*821E-08 

0. 3766-08 

17 

11474 

1362 

0.203E 

11 

0*l03E-07 

0.552E-08 

19 

1L474 

1372 

0*175E 

07 

0.247E-10 

0.5B7E-12 

20 

11474 

1372 

0.351E 

or 

0.5l7E-i0 

0.U7E-11 

21 

11474 

1372 

0.530E 

07 

0.752E-10 

0. 1416-11 

22 

11474 

1372 

0.123E 

08 

0* 168E-09 

0.2B2E-U 

23 

11474 

1372 

0,192E 

08 

0.261E-09 

0.4826-11 

24 

11474 

1372 

0.367E 

08 

0.458E-09 

0.799E-11 

25 

11474 

1372 

0.715E 

08 

0.679E-09 

0.132E-10 

26 

11474 

1372 

0^141E 

09 

0.869E-09 

0.200E-10 

27 

11474 

1372 

0.296E 

09 

0. 113F-08 

0i345E-10 

28 

11474 

1372 

0.567E 

09 

O*148E-0a 

0.674f^-l0 

29 

11474 

1372 

0.142E 

10 

0.261F-08 

0.2846-09 

30 

11474 

1372 

0-272E 

10 

0.409E-08 

O.034F-O9 

31 

11474 

1372 

0W35E 

10 

0.H0lE‘08 

0.341E-08 

32 

11474 

13 72 

0.U7E 

11 

0*1276-07 

0. 7056-08 

33 

11474 

1372 

0* 160E 

ll 

0*1696-07 

0. 108E-07 

34 

11474 

1372 

0.2030 

ll 

0.218E-07 

0.155E-07 

36 

11474 

1302 

0.175E 

07 

0.3346-10 

0.106E-11 

37 

11474 

1382 

0.351E 

or 

0.698E-10 

0.153E-11 

38 

11474 

1382 

0.530E 

07 

0.1026-09 

0.247E-11 



Table 

B-1. 

Turbid Water 

Data. 


PHHI MATTS) 


PHV(MATTS) 
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PARTICLE 

SAMPLE COWCENTRATION 


NUMBER 

DATE 

CODE 

(NUMBFR/MLl 

PVV(WATTS) 

PVH(M4TTS» 

PHH( WATTS) 

39 

11474 

1382 

0.881E 

07 

0.172E-09 

0.364E-11 


40 

11474 

1382 

0.123E 

08 

0.233E-09 

0.529E-H 


41 

11474 

1382 

0.192E 

08 

0.348F-09 

0.646E-U 


42 

11474 

1382 

0.367E 

08 

0.560E-09 

0 . ii 2 r-io 


43 

11474 

1382 

0.715E 

08 

O.0O7E-O9 

0.201E-10 


44 

11474 

1382 

0.141E 

09 

0.105E-08 

0.300F-10 


45 

11474 

1382 

0.296E 

09 

0. 138E-08 

0.593E-10 


46 

11474 

1382 

0.5676 

09 

0- 195E-08 

0- 138E-09 


47 

114 74 

1382 

0.142E 

10 

0.350E-08 

0.563E-09 


49 

11474 

13B2 

0.2T2E 

10 

0.636E-08 

O.214E-0B 


49 

L 1 4 74 

i38e 

0.7356 

10 

0.136F-07 

0. 799E-0B 


50 

11474 

1382 . 

0. 117E 

11 

0.202F-07 

0. 136F-07 


51 

11474 

1382 

0.160E 

ll 

0.272E-07 

0.204F-07 


52 

11474 

1382 

0.203E 

11 

0.324E-07 

0.258E-07 


54 

11574 

2372 

0.454E 

09 

0.634E-U 

0.235E-12 

0.752E-U 

55 

11574 

2372 

0.907E 

09 

0. 1416-10 

0.7056-12 

0.1416-10 

56 

11574 

2372 

0.136E 

10 

0.211E-10 

0.164E-U 

0.2116-10 

57 

11574 

2372 

0.227E 

10 

0.34BE-10 

0.352F-11 

O.340F-1O 

58 

11574 

2372 

0.31BE 

10 

0.4 ?0F-10 

0.493F-IL 

0.484E-10 

59 

11574 

2372 

0.499E 

10 

0.745E-10 

0.B22E-11 

0.750E-10 

60 

11574 

2372 

0.863E 

10 

0. 12BE-09 

0. 153E-10 

0. 1266-09 

61 

11574 

2372 

0.132E 

11 

0. 195F-09 

0-223E-10 

0.1956-09 

62 

11574 

2372 

0.222F 

11 

0.335E-09 

O.364E-10 

0.329E-09 

63 

11574 

2372 

0.3I3F 

a 

0.47ie-09 

0.505E-10 

0.4586-09 

64 

11574 

2372 

0.449E 

11 

0.6476-09 

0.74GF-10 

0.669E-09 

65 

11574 

2372 

0.722E 

11 

0. 104E-08 

0. U6E-09 

0. 1046-08 

66 

11574 

2372 

0.994E 

u 

0. 135E-08 

0. 154E-09 

0,1336-08 

67 

11574 

2372 

0.168E 

12 

0. 192E-08 

0.222F-09 

0.1916-08 

68 

11574 

2372 

0.304E 

12 

0.269F-08 

0.323E-09 

0.2606-08 

69 

11574 

2372 

0.576E 

12 

0.323E-08 

0.443E-09 

0.330E-08 

70 

11574 

2372 

0.112E 

13 

0.420E-08 

0-669E-09 

0.422E'-08 

71 

11574 

2372 

0.248E 

13 

0.563E-08 

0- 124E-08 

0.5636-08 

72 

11574 

2372 

0.657E 

13 

0.974E-08 

0.334E-08 

0. 9866-08 

73 

11574 

2372 

0. 107E 

14 

0.131E-07 

0.564E-08 

0.1316-07 


Table B-1, Turbid Water Data Continued, 


. PHV(WATTS) 


0*235£-l2 
0.470E-12 
0 , 7056-12 
0.21UE-U 
0. 31996-11 
0.799E-11 
0.1676-10 
0.26aE-l0 
0,^;a4E-l0 
U.7)07E-lO 
0.978E-10 
0. 1159E-09 
0.206E-09 
0.293E-09 
0.4226-09 
0.^B7E-09 
0.840E-09 
0.146E-08 
0.371F-08 
0.6116-08 
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SAMPLE 

NUMBER 


PARTICLE DYE 

CONCENTRATION CONCENTRATION 

(MG/L) 


DATE CODE 

11774 1372 

11774 1372 

11774 1372 

11774 1372 

11774 1372 

11774 1372 

U774 1372 

11774 1372 

11774 1372 


(NUMBER/ML) 

0*437E 10 
0.437E 10 
0.437E 10 
0.437E 10 
0.437F 10 
0.437E 10 
0.437E 10 
0.437E 10 
0.437E 10 


0. 

OOOE 

00 

0. 

809F 

00 

0. 

162E 

01 

0. 

4856 

01 

0. 

U3E 

02 

0. 

243E 

02 

0. 

502E 

02 

0. 

IIOE 

03 

0. 

2I4E 

03 


PVV( WATTS ) 

0.564E-08 
0. 536E-08 
0.517E-08 
0.482E-08 
0.421E-08 
0. 336E-08 
0.247E-08 
0.150E-08 
0.822E-09 


PVH( WATTS) 

0. 160E-08 
0. 132F-08 
0. 123E-08 
0. 105E-08 
0, 627E-09 
0.571E-09 
0. 343E-09 
0. 160E-09 
0.564E-10 


Table B-1. Turbid Water Data Continued. 
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PARTICLE OIL 


SAMPLE 

NUMBER 

DATE 

CODE 

CONCENTRATION 

(NUMBER/ML) 

OIL 

TYPE 

THICKNESS 

(MHI 

PVVIMATTSI 

PVHIHATTSI 

1 

U774 

1372 

Q,175E 

10 

NO 2 FUEL 

0.000 

0. 8936^08 

O.388E-O0 

2 

U774 

1372 

0.175E 

10 

N02 FUEL 

0.050 

0.728E^08 

0.2 94E-0:a 

3 

11774 

1372 

0. 175E 

10 

N02 FUEL 

O.lOO 

0.5 176-08 

0.2116-08 

h 

11774 

1372 

0.175E 

10 

N02 FUEL 

0.500 

0.172E-08 

O.601E-O9 

5 

11774 

1372 

0. 175E 

10 

N02 FUEL 

l.lOO 

0.3106-09 

0.1326-09 

6 

11774 

1372 

0. 175E 

10 

N02 FUEL 

2.500 

0.211E-10 

0.822E-IL 

7 

11774 

1372 

0.175E 

10 

SAE 30 

0.000 

0.8936-08 

0.399E-d8 

8 

11774 

1372 

0*175E 

10 

SAE 30 

0.050 

0.8S9E-08 

0.3766-08 

9 

11774 

1372 

0^175E 

10 

SAE 30 

0. 100 

0.8466-08 

0.352E-08 

10 

11774 

1372 

0*175E 

10 

SAE 30 

0.500 

0.8226-03 

0.3296-08 

11 

11774 

1372 

0.175E 

10 

SAE 30 

L. 100 

0.7756-08 

0. 3176-08 

12 

11774 

1372 

0.175E 

10 

SAE 30 

2.500 

0.728E-08 

0.294E-08 

13 

11774 

1372 

0*1756 

10 

SAE 30 

3.500 

0.6586-08 

0. 2706-08 

14 

11774 

1372 

0* 175E 

10 

SAE 30 

5.500 

0.5076-03 

0.235E-08 

15 

11774 

1372 

0.175E 

10 

GASOLINE 

0.000 

0*8466-08 

0.364E-08 

16 

11774 

1372 

0*175E 

10 

GASOLINE 

0.050 

0.0466-00 

0.3646-08 

17 

11774 

1372 

0.1T5E 

10 

GASOLINE 

O.lOO 

0.8466-00 

0.364E-08 

18 

11774 

1372 

0*175E 

10 

GASOLINE 

0.500 

0. 8226-08 

0.352E-08 

19 

11774 

1372 

0* 175E 

10 

GASOLINE 

1.100 

0.8116-08 

0.343E-08 

20 

11774 

1372 

0. 1756 

10 

GASOLINE 

2.500 

0. 7996-08 

0.329E-08 

21 

11774 

1372 

0*1756 

10 

GASOLINE 

3.500 

0.794E-08 

0.324E-08 

22 

11774 

1372 

0.1756 

10 

GASOLINE 

5.500 

0. 7876-08 

0.3206-08 

23 

11774 

1372 

0*1756 

10 

KEROSENE 

0.000 

0.8696-08 

0.376E-08 

24 

11774 

1372 

0.175E 

10 

KEROSENE 

0.050 

0.0466-08 

0.357E-08 

25 

11774 

1372 

0.1756 

10 

KEROSENE 

0.100 

0.8346-08 

0.3526-08 

26 

11774 

1372 

0.1756 

10 

KEROSENE 

0.500 

0.8226-08 

0.3506-08 

27 

U774 

1372 

0. 1756 

10 

KEROSENE 

I. 100 

0.8116-00 

0*341E'^08 

28 

11774 

1372 

0*1756 

10 

KEROSENE 

2.500 

0.8066-08 

0.329E-0B 

29 

11774 

1372 

0.1756 

10 

KEROSENE 

3.500 

0.804E-08 

0.3276-08 

30 

11774 

1372 

0. 175E 

10 

KEROSENE 

5.500 

0. 7076-08 

0.3I7E-^08 

31 

11774 

1372 

0.397E 

u 

L CRUDE 

0.000 

0.376E-07 

0.305E^OT 

32 

11774 

1372 

0.397E 

11 

L CRUDE 

0.050 

0.9056-08 

0.T05E-08 

33 

11774 

1372 

0*3976 

ll 

L CRUDE 

O.lOO 

0.329E-08 

O.250E-O8 

34 

11774 

1372 

0.397E 

u 

L CRUDE 

0.200 

0.2476-09 

0. lQflE-09 

35 

11774 

1372 

0*3976 

11 

L CRUDE 

0.300 

0.1766-10 

0.1366-10 


Table B-2. Oil on Water Data. 
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SAMPLE 

NUMBER 

DATE 

CODE 

PVV(WATTS) 

PVHXWATTS) 

4 

7L473 

3X41 

0-493E-10 

0.106E-10 

14 

71473 

3131 

0.388E-10 

O-7526-ll 

23 

71473 

3121 

0.258E-10 

0.3526-n 

36 

71473 

3111 

0.223E-10 

0.352E-11 

54 

71573 

31U 

0,458E-10 

0.940F-11 

59 

71573 

3111 

0.435E-10 

0.7286-U 

60 

71573 

3U1 

0.4UE-10 

0.6llE-ll 

az 

71573 

3111 

0-540E-10 

0.113E-10 

345 

92573 

3251 

O,180E-O9 

0.294E-10 

2 30 

92573 

3251 

0. 153E-09 

0.141E-10 

326 

92573 

3251 

0,U7E-09 

O.IAIE^IO 

134 

92573 

3251 

0, U7E-09 

O.B22E-U 

112 

92573 

3251 

0*117E-09 

0.940E-1 1 

362 

92573 

3251 

□•022E-1O 

0,9406-U 

113 

92573 

3251 

0.9636-10 

0,IQ6E-10 

150 

92573 

3251 

0. 1126-09 

0.106E-10 

34 

92573 

3251 

0.822E-10 

0. 8226-11 

235 

92573 

3251 

0.681E-10 

0.940E-LI 

99 

LZ0173 

3252 

0.235E-09 

0. 38BE-10 

11 

120173 

3252 

0.329E-09 

0.705E-10 

57 

120173 

3252 

0.399E-09 

0.999E-10 

225 

120173 

3252 

0.587E-09 

0.200E-09 

24 

120173 

3252 

O.200E-09 

0.376E-10 

59 

120173 

3252 

0.282E-09 

0.5T6E-10 

41 

120173 

3252 

0,67QE-09 

0.25BE-09 

264 

12CU73 

3252 

0.7526-09 

0.305E-09 

217 

120173 

32 52 

0. 188E-09 

0. 3526-10 

165 

120173 

32 52 

0.846E-09 

0. 364E-09 

7fl 

120173 

32 52 

0. 70SE-09 

0.294E-09 

250 

120173 

3252 

0.634E-09 

0.223E-09 

90 

120173 

3252 

0.799F-09 

0. 329E-09 

134 

120173 

3252 

0. 1016-08 

0.5176-09 

56 

120173 

3252 

0.176E-09 

0.235E-10 

23 

120L73 

3252 

0.294E-09 

0. 5646-10 

5Z 

120X73 

3252 

0. 1366-09 

0.1416-10 

360 

120173 

3252 

0.940E-10 

0.1136-10 

352 

120173 

3252 

0*799E-10 

0.9406-11 

35 

120173 

32 52 

0.98TE-10 

0. 1136-10 

135 

120173 

3252 

0.164E-09 

O.2356-10 

70 

120173 

3252 

0. 2356-09 

0. 3996-10 

386 

120173 

3252 

0*2006-09 

0.329E-10 


PHH< WATTS) 

0*493E-10 

O.3a8E-’l0 

0*^47E'10 

0*2UF-10 

0.4466-10 


PHVCWATTS) 

0,ll7S-lO 

0.775E-11 

0.35^£-ll 

0^352e-ll 

0*940E-ll 


SUSPENOEO 


ToSBioiry 

TRANSMITTANCE 

SOLIDS 

<FTU> 

{%) 

IMG/L) 

7.0 

93.7 


8.0 

93.5 


8.0 

92-5 


9.0 

94.0 

10. 0 

28.0 

75.0 

J8.0 

21.0 

81. 0 

31. 0 

22. 0 

81.0 


32.0 

73.0 

49.0 

26.0 

84.0 


22.0 

86.0 


20.0 

87.0 


18.0 

88.0 


19.0 

67.0 

— 

18.0 

69.0 


Id.o 

38.0 


IT.O 

89.0 


16.0 

90.0 


15.0 

90.0 


46.0 

61.0 

J 13.0 

82.0 

63.0 

141.0 

215.0 

40.0 

266.0 

350.0 

19.0 


42.0 

61.0 

125.0 

75.0 

60.0 


450.0 

10.5 

603.0 

575.0 

6.5 

814.0 

25.0 

66.5 



650.0 

5.0 


550.0 

7.2 


500.0 

10. □ 

430.0 

625.0 

4.5 


700-0 

2.8 

2056.0 

28.0 

84.0 


69.0 

70.5 


16.0 

89.0 

— — 

12.0 

95-0 


li.O 

96.0 

62.0 

12«0 

96.0 

— — 

27.0 

89.5 

63.0 

44.0 

77.5 


45.0 

81.5 




Table B-3. Field Data from Natural Waterways 



